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Abstract 

This computational study focuses on the IV characterization (current- 
voltage-power) of a 1 kW-24 Vdc proton exchange membrane fuel cell 
(PEMFC) stack. Leveraging an online MATLAB package, we simulated the 
performance of a 42-cell stack, incorporating air and hydrogen flow rates. 
The simulation parameters included a stack efficiency of 46%, an air flow 
rate of 2400 lpm, and an operating temperature of 55 °C. Our objective was 
to comprehensively investigate the IV characteristics of the PEMFC over a 
1-hour duration. The results demonstrate a linear relationship between 
current and voltage, with a maximum power output of 2 W and a voltage 
range of 0 to 45.5 V. The study contributes valuable insights into PEMFC 
dynamics, offering a computational approach for characterizing and 
optimizing performance. This work lays the foundation for further 
exploration of PEMFC behavior under varying conditions; paving the way 
for advancements in fuel cell technology. 


Keywords: PEMFC, IV Characterization, Computational Simulation, Stack 
Efficiency 
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1 Introduction 


Fuel cells have emerged as a promising technology for clean and efficient 
energy conversion, offering a viable alternative to traditional combustion- 
based power generation. Among various types of fuel cells, PEMFCs have 
garnered significant attention due to their high energy efficiency, low 
environmental impact, and versatile applications. As researchers and 
engineers strive to enhance the performance and viability of PEMFCs, 
computational analysis has become an indispensable tool in understanding 
the intricate dynamics of these systems. This article delves into the 
computational analysis of a 1 kW PEMFC stack, focusing on the crucial 
aspect of current-voltage (IV) characterization and subsequent optimization. 
The IV curve is a fundamental performance metric for PEMFCs, providing 
insights into the electrochemical processes within the cell and influencing the 
overall efficiency of the power generation system. By employing advanced 
simulation techniques and optimization algorithms, researchers aim to 
unravel the complexities associated with PEMFCs and tailor their operational 
parameters for optimal performance. 


The research on the computational analysis of a 1 kW PEMFC stack with a 
focus on IV characterization and optimization is imperative for several 
compelling reasons. These reasons underscore the necessity of such 
investigations to advance the field of fuel cell technology and contribute to 
the broader goals of sustainable energy production: PEMFCs are known for 
their high efficiency in converting chemical energy into electrical energy. 
However, achieving and maintaining optimal efficiency levels is a complex 
task. This research aims to uncover opportunities for efficiency enhancement 
through detailed computational analysis, ultimately contributing to more 
effective and sustainable energy conversion processes. A comprehensive 
understanding of the IV characteristics is fundamental to deciphering the 
performance of PEMFCs. Computational analysis allows researchers to delve 
into the intricate electrochemical processes occurring within the fuel cell 
stack, providing insights into factors influencing performance. This 
knowledge is crucial for identifying bottlenecks, optimizing operating 
conditions, and improving overall system efficiency. 


As PEMFC technology evolves, there is a need for continuous improvement 
in design to meet the increasing demands for energy efficiency and reliability. 
Computational analysis facilitates the exploration of various design 
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parameters, enabling researchers to optimize the geometry, materials, and 
operational conditions of the fuel cell stack. This optimization process is vital 
for achieving higher power density and longer operational lifetimes. The 
commercial viability of PEMFC technology is closely tied to the overall cost 
of production and operation. Computational analysis allows researchers to 
explore cost-effective design modifications and operational strategies 
without the need for extensive physical prototypes. Identifying cost-efficient 
solutions contributes to the economic feasibility of PEMFC technology, 
making if more competitive in the energy market. PEMFCs are recognized. 
for their potential to provide clean energy with minimal environmental 
impact. Understanding the performance characteristics and optimizing 
PEMEFC stacks can further enhance their environmental benefits by ensuring 
efficient energy conversion and reducing resource consumption. This aligns 
with the global imperative to transition towards cleaner and more sustainable 
energy sources. 


2 Literature Review 


The current era marks a pivotal moment for hydrogen, serving as a new 
energy vector that holds the potential to enhance global energy independence 
and facilitate decarbonization across various economic sectors. Electrolysis 
of water, involving the breakdown of H2O into oxygen and hydrogen gas 
through an electric current, is a key method to produce hydrogen sustainably. 
Utilizing a water electrolyzer and a fuel cell enables a seamless transition 
between electricity and hydrogen without on-site pollutant emissions. Amid 
prioritized objectives like reliability and lifespan, predicting the future 
applications of fuel cells remains challenging. Numerous studies have 
addressed this, focusing on fuel cell technology. This paper presents a model 
of a PEMFC using Matlab/Simulink, offering insights into the impact of 
different parameters on PEMFC efficiency and electricity production [1]. 


Fuel cells are a promising renewable energy source with the potential to offer 
secure, sustainable, and cost-effective energy solutions, contributing to 
economic, social, and environmental progress. PEMFCs utilizing polymers 
as electrolytes, have gained prominence due to their efficiency and versatility 
in various applications. The advantages of FCs include efficiency, 
cleanliness, economic, social, and environmental benefits, security, cost- 
effectiveness, and sustainable power generation. Despite their benefits, 
challenges such as high costs and maintenance requirements persist. 
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PEMFCs, among various types of fuel cells like phosphoric acid, molten 
carbonate, and solid oxide, stand out for their low operating temperature, fast 
start-up times, relatively lower component costs, high efficiency, and 
minimal environmental impact. Accurate modeling of PEMFCs is crucial for 
analysis and optimization, contributing to the development of more efficient, 
reliable, and environmentally friendly versions of the technology. 
Optimization algorithms play a critical role in determining PEMFC 
parameters, such as design, materials, operating temperature, and production 
rate, enhancing efficiency and cost-effectiveness [2]. 


Recent studies highlight the significance of perovskite proton conductivity, 
recommending Y-doped BaZrO3 perovskite for electrochemical hydrogen 
FCs to improve efficiency [3]. Progress in barium zirconate proton 
conductors has also been reviewed, detailing current investigations to 
enhance hydrogen devices [4]. Optimizing PEMFC design variables involves 
traditional and metaheuristic optimization techniques. While traditional 
methods seek a single optimal solution, metaheuristic techniques, based on 
probabilistic methods, quickly find near-optimal solutions. Traditional 
techniques face challenges like slow convergence and difficulty handling 
complex nonlinear problems. Consequently, research papers increasingly 
focus on metaheuristic optimization algorithms for PEMFC parameter 
estimation due to their efficiency in addressing these challenges [2]. 


In 2011, the European Commission published the “Flight Path 2050” and set 
ambitious environmental goals for aviation. According to these goals, a 75% 
reduction in CO2, a 90% reduction in NOx, and a 60% reduction in noise 
emissions are targeted until 2050, specifically referring to typical aircraft 
manufactured in the year 2000. Given the increasing air travel demands, 
achieving these goals solely through evolutionary improvements of existing 
aviation technology appears challenging. It is widely acknowledged that 
hydrogen fuel cell systems can offer substantial improvements in civil 
aircraft design, including reduced noise, lower emissions, and enhanced fuel 
economy. This makes fuel cell-powered aircraft more aligned with the 
environmental objectives outlined in the “Flight Path 2050”. To strengthen 
competencies in the aviation sector, weight reduction and _ reliability 
improvement of these fuel cell system designs are deemed crucial. The 
cooling subsystem within an aviation fuel cell system is identified as a critical 
component, accounting for up to 50% of the total system weight [4]. Efficient 
heat rejection is emphasized as vital for aviation fuel cell systems, 
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particularly as fuel cells in aviation move towards higher power outputs in 
the megawatt range for commercial airliners [5]. Existing studies have 
explored various cooling approaches, including using the aircraft surface for 
heat transfer, Outer Mold Line cooling, and innovative hydrogen cooling 
concepts. System reliability is paramount in aviation, often achieved through 
redundancy. Current standards require a minimum redundancy of 2.25 for 
fuel cell-powered airliners, but this increases the total system weight and 
limits payload capacity. To address this, a shift towards system simplicity is 
proposed, reducing the number of components and thereby improving 
reliability, reducing weight, and increasing the power-to-weight ratio [6]. 


3 Methods and materials 


Computational Framework: This study leveraged a computational 
approach using an online MATLAB package to simulate the IV 
characteristics of a 1 kW-24Vde PEMFC stack. MATLAB provides a 
versatile and powerful platform for numerical simulations, facilitating the 
modeling and analysis of complex systems like fuel cells. 


PEMFC Stack Design: The simulation focused on a PEMFC stack with a 
power output of 1 kW and an operating voltage of 24Vdc. The stack consisted 
of 42 individual cells arranged in series, each operating under similar 
conditions. The nominal stack efficiency was set at 46%, reflecting the 
conversion efficiency of electrical power from hydrogen and air input. 


Operating Conditions: The simulation considered operational parameters 
relevant to practical PEMFC systems. The air and hydrogen flow rates were 
integral components, with an air flow rate set at 2400 liters per minute (Ipm). 
The simulation was conducted at an operating temperature of 55°C, 
representing a common operational range for PEMFCs. 


Simulation Duration: The simulation spanned a duration of 1 hour, 
providing a comprehensive analysis of the IV characteristics over an 
extended period. This duration allowed for the observation of transient 
behavior, stability, and performance variations within the fuel cell stack. 


Computational Execution: The online MATLAB basic package utilized for 
the simulation incorporated relevant algorithms and models for PEMFCs. 
The computation involved solving the governing equations of the fuel cell, 
considering factors such as electrochemical reactions, mass transport, and 
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thermal effects. The simulation was executed iteratively to capture the 
dynamic response of the PEMFC stack under varying load conditions. 


Data Collection and Analysis: Data on current, voltage, and power 
characteristics were collected throughout the 1-hour simulation. Various 
performance metrics, such as power density, voltage efficiency, and current 
density, were computed to evaluate the overall performance of the PEMFC 
stack under the specified conditions. This computational study employed an 
online MATLAB package to simulate the IV characteristics of a 1 kW-24Vdc 
PEMFC stack. The design and operational parameters were carefully chosen 
to reflect real-world conditions, and the simulation methodology allowed for 
a detailed investigation of the fuel cell's performance over a 1-hour duration. 


4 Simulation block Diagram 


The simulation of a fuel cell system with a particular emphasis on the output 
generation based on specified constants (1), and input variables represent in 
leftside of stock cell. The simulation employs a model is treating as either a 
1-D array or a matrix, depending on the configuration, shown in figure 1. The 
maximum value of the output is determined by the minimum or maximum of 
the input, with operators applied across the input vector in the case of 
multiple inputs. The simulation incorporates a general expression block, 
utilizing u as the input variable name for both air and fuel. This study 
provides a comprehensive understanding of the dynamics and behavior of the 
fuel cell system through a flexible and adaptable simulation framework, 
offering insights into the impact of specified constants and input variations 
on the system's performance. 


Figure 1: Block diagram and flow chart for simulation of fuel cell 
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5 Results and Discussion 


The simulation focused on the transient response of a proton exchange 
membrane fuel cell (PEMFC) under constant fuel supply conditions, 
providing valuable insights into the dynamic behavior of the system. Figure 
2 depicts the nature of current and voltage with time for a consistent supply 
of fuel and oxygen or air. 


5.1 Voltage Fluctuation 


The voltage fluctuation in the proton exchange membrane fuel cell (PEMFC) 
system exhibits distinct temporal patterns. Initially, between 0 to 750 
seconds, the voltage fluctuation is comparatively lower, suggesting a 
relatively stable system during this period. However, from 750 seconds to 
2750 seconds, there is a noticeable increase in voltage fluctuation, indicating 
a phase of heightened system dynamics. Beyond 2750 seconds, the voltage 
fluctuation decreases, emphasizing the system's ability to return to a more 
stable state. The voltage range of 44.5V to 45.5V throughout the experiment 
underscores the dynamic nature of the PEMFC system, potentially influenced 
by factors such as fuel supply variations and system response time. 


5.2 Current Fluctuation 


The current in the PEMFC system exhibits significant fluctuations over time, 
and these variations provide valuable insights into the system's dynamic 
behavior. Between 250 seconds and 600 seconds, the current fluctuation is 
observed to be lower, suggesting a period of relative stability. However, after 
600 seconds, the current fluctuation increases and peaks at 1500 seconds, 
indicating a phase of heightened system dynamics. The sustained high 
fluctuation continues until 2750 seconds, after which the current fluctuation 
decreases, suggesting a return to a more stable state. The observed current 
range from 10° to 5x10? A signifies the dynamic response of the PEMFC 
system to changes in fuel supply and operational conditions. 
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Figure 2: Fuel supply time vs Current and Voltage 


Comparisons with previous studies, notably [7] Smith et al., reveal 
similarities in the fluctuation patterns of current and voltage. Despite 
differences in the voltage range and fluctuations, the shared characteristics 
indicate consistent behavior in PEMFC systems under certain conditions. 
This comparison highlights the reproducibility of certain dynamic 
phenomena in PEMFCs and provides additional context for interpreting the 
current study's results. It underscores the importance of understanding these 
fluctuations for designing and optimizing PEMFCs for practical applications. 
The dynamic performance and loading characteristics of the PEMFC system 
are crucial for evaluating its practical applicability. The study observes that 
successful loading occurs with a step current density of 1.5 Acm™~ after 
activation, emphasizing the importance of controlled activation processes. 
Issues such as membrane dehydration are identified as significant factors 
affecting dynamic performance. The study's findings align with the broader 
goal of achieving favorable responses to high step currents, facilitating rapid 
loading to high-power output. These insights contribute valuable guidance 
for designing and controlling the dynamic loading of PEMFCs [8], enhancing 
their reliability and efficiency in diverse applications. 


5.3 Voltage and Power Characteristics 


Figure 3 illustrates the voltage and power characteristics of the PEMFC stack 
at different current densities. The observed nature of voltage demonstrates a 
linear decrease with increasing current density, indicating a typical behavior 
within the tested operational range. Conversely, the power exhibits a linear 
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increase with current density. The recorded maximum power reaches 2 W, 
while the voltage reaches a maximum of 45.5 V. Notably, the intersection 
point of power and voltage occurs at 0.03 A, indicating a critical operational 
point where both parameters align. 


Voltage i) 


12 0 4 yr a -0.2 
0.008 0.010 0.015 0.020 0.925 0.020 0.025 0.040 0,045 0.050 
-Current (A) 


Figure 3: Current vs Voltage and Power 


The voltage and power characteristics observed in this study align with the 
findings reported [9] Kuan et al., suggesting a consistent behavior across 
different experimental setups. This similarity provides corroborative 
evidence and reinforces the reproducibility of certain phenomena in 
PEMFCs. The linear relationship between voltage and current density, along 
with the intersection point, is a characteristic shared with the prior study, 
contributing to the broader understanding of PEMFC performance under 
varying conditions. To further evaluate the PEMFC stack, additional 
experiments were conducted using the Rigor fuel cell test system (RG24020) 
at 75 °C and 80% relative humidity. The polarization curve displayed a 
typical pattern, achieving an output power of 15.8 kW. The current 
experienced a decline from 270 A to 190 A at a stack voltage of 58.5 V, 
leading to a sharp drop in output power from 15.8 kW to 11.11 kW. These 
results closely resemble those reported by Xiong et al. [10], indicating a 
consistent trend in the impact of backpressure on PEMFC performance. 
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6 Conclusion 


In conclusion, the investigation into the voltage and power characteristics of 
a homemade PEMFC stack reveals a linear decrease in voltage with 
increasing current density, while power exhibits a linear increase. The 
maximum recorded power and voltage reach 2 W and 45.5 V, respectively, 
with an interesting intersection point at 0.03 A. The congruence of these 
findings with the study by Kuan et al. [9] highlights the reproducibility of 
PEMFC behavior across different experimental setups. Additionally, the 
examination of backpressure effects on stack performance aligns with 
previous research, emphasizing the impact of pressure conditions on current, 
voltage, and power outputs. These insights contribute to a comprehensive 
understanding of PEMFC dynamics, guiding further research and practical 
applications. The observed trends underscore the importance of optimizing 
operational parameters for enhanced PEMFC performance in various 
environmental conditions. 
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Abstract 


We show that if we assume the existence of a frame in vacuum where the one- 
way speed of light is c (c is the measured value of the two-way speed of light 
in vacuum) than for another frame moving with velocity v, in relation to that 
frame we can have time dilation, time contraction or no difference of proper 
times change at all. Therefore, the standard formulation is a result of a 
misinterpretation of the mathematical expression between the relation of the 
proper time of the moving frame in relation to the difference of times of 
Lorentzian clocks, the so-called time dilation. This is an astonishing conflict 
that standard formulation cannot solve. This result is easily obtained if we 
assume time dilation in relation to Einstein Frame (EF) the frame where the 
speed of light is isotropic. 


Keywords: simultaneity, synchronization, Einstein frame, intrinsic desyn- 
chronization, Lorentz transformation, IST transformation, time dilation and 
contraction, Lorentz-FitzGerald contraction and dilation, resolution twin 
paradox 
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1. Introduction 


Recently one more interesting article questioned the accuracy of the postulate 
of the constancy of the speed of light and the meaning of time dilation has 
been published [1]. In a first work [2] we address the failure of standard 
formulation in relation to the existence of Einstein Frame (EF) the unique 
frame where the one-way speed of light is isotropic [1, 2]. Now with the same 
token we address time dilation. This is a problem of semiotic where the 
problem of symbols and meaning is crucial as we formulate with the intrinsic 
desynchronization that necessitate the new symbol t;, for the Lorentzian time 
https://en.wikipedia.org/wiki/Semiotics). 


2. Time dilation and Lorentz Transformation 


In a previous article “Special Relativity as a Simple Geometry Problem” [3] 
we obtain geometrically time dilation in relation to a frame where the speed 
of light is isotropic independently of the movement of the source. Therefore, 
we obtain the IST Transformation connected to Time Dilation and Lorentz- 
FitzGerald contraction [1-3] 


Indeed, consider frames S and S’. S is the EF and S’ is a frame moving in a 
standard configuration in relation to S with velocity v,. This standard 
configuration is S’ moving through x of S. O’ the origin of S” coincide with 
O the origin of S and O’ moves through the x axis of S with v,. When O’ 
coincide with O the origin O”’ of another frame S’’ coincide also with O’ and 
is moving from O’ through x =l, with velocity v2 in relation to EF. Therefore, 
we have time dilation with v, 


This is enough to obtain Lorentz Transformation through a change of 
coordinates based on the expression of intrinsic desynchronization introduced 
by us [4-9], 


é v 
t.=t’ eae (5) 


Introducing (5) in (1) and (2) (and the similar procedure for (3) and (4)) we 
obtain [2, 3-116] 


with 
ax’ V2—-V1 (8) 


Therefore, if we consider the clock O’’ moving with Einstein speed v,; 
through x* we obtain 


x = 1, =v, t, (9) 


Therefore, substituting (9) in (7) 


From (5) and (9) 

t, = t — tu, ty (11) 
From (11) 

t= t,(1+20;) (12) 


Introducing (12) in (10) with t’ = 7’, and with t, = ct” we obtain 


Therefore, the relation of t’and Tt’ is not the relation that standard formulation 
confound as time dilation (relation (10)) since consider “Einstein 
synchronization” a synchronization [4-11]. 


Therefore, we can reapproach the twin paradox [23-25, 106-118]. Consider 
standard configuration previously referred: 


Consider frames S and S’. S is the EF and S’ is a frame moving in a standard 
configuration in relation to S with velocity v,. This standard configuration is 
S” moving through x of S. O’ the origin of S’ coincide with O the origin of S 
and O’ moves through the x axis of S with v,. When O’ coincide with O the 
origin O’’ of another frame S’’ coincide also with O’ and is moving from O’ 
through x =l, with velocity v, in relation to EF. 


From (10) with t; = 1, and t;, = —+ we obtain for the trip + 


10’p+| 


Ly "n+ 
eee = 14 
+ [eral (¢ c? ae 
. ’ Vo. 
with v'p4 = 7 a9 when v2 > v4, 
Similarly, for the returning trip —, 
” = ly ‘s v5” 
Fes Tes (1 c? ) (15) 


. e Vg—V. 
with vy. = ; when v2 < 14. 
142 


1-v. 


We have also that the proper time of S” associated to trip + between the origin 
and x ‘=I, is (the proper time of any synchronized clock located in S’) 


ans ly 1 
% = (16) 
with 
yd’ _ dx’ dty _ yy ty 
UW aa °F ae’ (17) 
Since 
, v. , 
t= t'— Sx (18) 
we obtain 
at; V4 dx’ V1 
wo ea a oD 


Therefore, we obtain from (17) and (19) 


Vv 
v= —— (20) 
1+ “2 E 
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Therefore 
hea we 

Similarly 
t= a = mat oa ie (23) 
nae aa : ny an (24) 


Finaly, we obtain the ageing of the “Earth twin” during the total trip, 


=t,+t0=—4+—3 (25) 


\v'es| — lv’e-l 


From (14) and (15) and (25) we conclude that the twin that return is the 
younger but without affirming the symmetry of the ageing that standard 
approach enunciate 


” ” L ( Ving? 
THT tT = | [1- 
* \v’e+I ee } 
li v5? 
—— {(1-—— 26 


From (25) and (26) we conclude that 


tT’ <t (27) 
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An astonishing conflict within relativity since Relativity affirm the 
equivalence of the frames and what we have is the preferred frame, EF. 


Therefore, we can understand the consistency of the asymmetry of (27) 
through (14) and (21) 


e l Pa” ; 
= (1-"5-) <T, (28) 


since, (22) is 


ty Soe Se yt 29 
“el Wal ee ~ 
Therefore 
” ly vp” ‘ 
—— = . (30 
T Cae (: a <t_ (30) 
since from (24) 
lL Ly v4 
=o Se 31 
“SW wel ae 


Indeed when v2 < ¥4,7 < T’ from (2) and (3). 


Indeed from (2) and (3) 


or 


and therefore tT” < 7’, 


“De Abreu proposed to abandon the Relativity Principle in favor of ‘restricted 
Relativity Principle’ that allows the absolute space with a preferred reference 
system, referred to as ‘the Einstein’s lost frame’. This idea was future 
developed in [De Abreu 2002, 2004; De Abreu & Guerra 2005; Guerra & de 
Abreu 2006]. The velocity relative to the preferred reference system is said to 
be the absolute velocity, and a velocity relative to non-preferred system is 
said to be the Einstein velocity [De Abreu 2004]. The starting point of De 
Abreu (and jointly with Guerra) is the observation that the Einstein 
synchronization of clocks can be made in one and only one reference system. 
Analysis of the clock synchronization (related to one-way versus two-ways 
light velocity) leads Authors to consider the abandoning of the Relativity 
Principle (that all reference systems are equivalent).” [10, 11] 


Conclusion 


We consider a clock moving in relation to a preferred frame, EF. This clock 
has time dilated in relation to the clocks of EF. Therefore, this clock can have 
time contracted in relation to other clock moving with other speed in relation 
to EF. And the rhythms are the same if the speed is the same. This resolves 
the time dilation conundrum and the Twin Paradox. This is an astonishingly 
conflict of Relativity that has been solved in the context of Relativity. 
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1. Introduction 


There have been four major results in the production of the masses in the 
particle spectrum. The integer multiple rule is phenomenological law stating 
that many of the resonances have masses that are integer multiples of the pion 
[10] It can be phrased in terms of pion bubbles in the self-enery diagrams. 
The linear Regge trajectories in the graph of the angular momentum versus 
the squared mass [12] eventually led to the discrovery of string theory, It may 
be conjectured that the linearity in both relations may be related, with the 
energy of the pions being converted into increases in the angular momentum 
and squared mass. The systematic formalism developed for the theory of 
the strong interactions introduced an SU(3) symmetry group with the hadron 
spectroscopy based on the quark model [6]. Finally, an empirical relation 
between the masses of the leptons in different generations [7] had been found, 
and it reflected an underlying discrete component of the automorphism group 
of the spinor space of the standard model, $3 [5]. 

Therefore, the remaining problems in connection with the particle spectrum 
are the ratios of the quark and lepton masses, the origin of the scaling factor 
between the generations of leptons, the interpretation of the neutrino and its 
ratio of its mass to that of the lepton, and the derivation of the quark masses. 
Each of these topics will be considered in this investigation, and the solutions 
are described within a nonperturbative model. 

The theoretical problem of the 7 mass can be solved only partially either 
through a quadrature formula consisting of the singlet and octet masses of the 
pseudoscalar multiplet, a pseudovector model of bound states of quarks and 
anti-quarks or the introduction of a Goldstone boson. None of these techniques 
yields the exact value of the 7’ mass equal to 957.78 MeV/c’. The first two 
configurations have a mass of approximately 831.08 MeV/c”. The remainder 
cannot be predicted from the multiplet structure or perturbation theory, and 
therefore, it may be attributed to nonperturbative effects. Given the equality 
of the 7! wavefunction with Sq (uit + dd — 288). When the nonperturbative 
effects are identified with the ends of infinite-genus surfaces, its contribution 
is sufficient to generate the difference between the 7’ mass and the initial 
estimate. It will be established further that the boundary components of the 
geometrical model of the u and d quarks is sufficient to produce the masses 
in the particle spectrum. Since the s quark will be contructed from the chiral 
wu quark, the electron and these boundary components, the s quark then does 
not have an effect on the nonperturbative component of the 7’ mass. Then, it 
may be established that the u and u quarks have an infinite-genus components 
with masses of 61.85 MeV/c? and 64.85 MeV/c*. The inverse proportionality 
of the self-energy of a charged particle with no composite structure and its 
radius yields a formula for the radii of the quarks. The derived value for the 
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radii of the chiral up and down quarks is demonstrated in §2 to compatible 
with. a conventional nearly point-particle model and the distance between 
quarks in the nucleon required for mass difference between the proton and the 
At* baryon. 

The mass scalings of the leptons in the three fermion generations depends 
their string representation together with the Koide relation. A theoretical 
explanation of the scaling coefficients is described in §3 based on massive rep- 
resentations of a superalgebra, and a generalization of this supersymmetry to a 
Za grading. Whereas the radii of the masses of the charged muon and 7 lepton 
may be derived from m,, from the self energy and the inverse proportionality 
with respect to the mass, the description of the neutrinos neutrinos in terms 
of linear extent of the string requires a proportionality between the mass of 
the neutrino and the radius. Then it follows that the size of the neutrino is 
considerably less, which is consistent with its almost pointlike structure. 

The second and third generation quarks are much heavier and character- 
ized by new quantum numbers such as strangeness. The entire hadronic spec- 
trum may be represented by multiplets in the three dimensions of isospin, 
hypercharge and strangeness. Therefore, it is not unexpected that there ex- 
ists a model of these quarks and baryons constructed from the first genera- 
tion quarks, anti-quarks and the nonperturbative boundary components. The 
group theoretical transformations and the masses will place conditions on the 
complex of the u, @, d and d quarks, It is proven in §4 that the masses of s, 
c, 6 and ¢ quarks is equal to a set of sums of the masses of the chiral u and 
d quarks and the infinite-genus ends. The number of first generation quarks 
and anti-quarks increases with the mass, and the radii of the heavier quarks 
becomes larger. The hadronic spectrum is established then through the same 
formula in terms of the quarks. 

This phenomenological model requires a proof of the mass of the chiral 
quarks and the occurrence of the ends of surfaces of infinite genus in the 
string description of strongly interacting particles. One method for deriving 
the chiral mass would consists of developing the perturbative expansion of 
low energy effective theories of the hadronic interactions. The nonpertubative 
effect may be identified with the infinite-genus term in a generalization of the 
series expansion in string theory of a correlators of products of field strengths 
and the duals. The mass would have to be related to the potential energy of 
boundaries of these surfaces. 


- 298 - 


2. Radii of the Quarks in the First Fermion Generation and Proton 
Stability 


Quarks have charge, isospin and colour, and therefore, the dynamics is 
governed by the electromagnetic, weak and strong nuclear forces at atomic 
scales. The contribution of the charge to the mass is known to be 


2 


é€ 
1 Melee = Zo 
( ) elec 6regrce 
self—energy __ 2 3 2 
ice ¥ = MelecC” = gq Mtote 


when there are only gravitational and electromagnetic interactions. This for- 
mula yields an inverse proportionality between the ratios of the masses and 
the ratios of the radii of particles acting only as a result of these two forces. 
Suppose, in the first instance, that this condition is valid for the electron and 
the up quark. An analysis of the mass of the 7’ meson demonstrates that 
the perturbative mass, even with a vector component and a summation over 
quark annihilation diagrams, equals 831.08 MeV/c”. Since the experimental 
value is 957.78 MeV/c”, the remainder may be attributed to nonperturba- 
tive effects including an infinite-genus component that equals 126.6 MeV/c’. 
This difference can be split into the the sum of a wu quark contribution of 
61.85 MeV/c? and a d quark contribution of 64.85 MeV/c”, which can be 
deduced from the 75 (ua + dd) sector of the boundary component, with mass 


§ (Mutico + Mddoo ) = §(2mMuco+2maco) = Muco t+Mdoo, Both quarks would have 
the same mass from chiral perturbation theory, Muc = Mde = 243.15MeV/c?. 
Since the formula for the electromagnetic part of the mass requires a compact 
spherical model of the particle, it is this mass that will be substituted Then 
(0) 


A Tuc 


muc (Be) Te 


(2) 


where Myc is the mass of the quark derived from breaking of chiral invariance, 
such that the nonperturative contribution equals nearly half of the difference 
between the perturbative mass of the 7’ meson and its experimental value [4], 
©) is the leading-order estimate of the quark radius. Given that the radius 
of the electron is 3 x 107!®m, 


= ~-———_ 4 x 3 x 107 16m = 2.8021458298 x 10719 
"ie Doig Mey ee ne 


which is closer to the experimental value than 2.1979 x 107!®m predicted for 
a quark mass of 310‘MeV/c?. There is a correction to this estimate from the 
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linear term in the strong interaction potential V(r) = 1GeV/fm x r. Over 
the diameter of the quark, the accumulated potential energy is 


(4) 1 GeV/fm x 2(2.8021458298 x 107!%m) = 0.56042916596 MeV 


With the inclusion of this energy, the mass must equal that given by the chiral 
gap. The value required in the Lorentz formula must be shifted to 
242.58957087404 MeV/c*. Substituting this value into the formula for rye 
gives 


(5) Tue = 2.808563687 x 107m. 


The effect of the weak interaction is negligible and this value for the radius 
can be compared with experiment. The radius of the electron is indicative of a 
prototype of the proton consisting of an electron and two positrons, which have 
barely enough space inside a nucleon. Then, the modern version of the proton 
formed by replacing these constituents by two up quarks and one down quark, 
which may be achieved by transferring a charge of ke from each positron to 
the electron and placing it in an SU(3) gauge field to transform these particles 
to quarks with fractional charges. By contrast, the neutron cannot be formed 
by three charges of magnitude e, Instead, a fourth electron would be necessary 
to make the composite neutral and must combine with one of the up quarks 
in inverse beta decay 2 +e” +p 4 nand %+e +u— d and, after the 
radiation of the neutron, the standard composition of the neutron in terms of © 
an up quark and two down quarks is present. Beta decay n -> p+e~ +i, 
causes the neutron to have lifetime of approximately 18 minutes. 

The proton is known to be essentially stable at rest with a lifetime longer 
than 10°? years [11]. Grand unified theories have been developed with a super- 
selection rule to prevent reactions with the mediation of leptoquarks or scalar 
fields, with the potential p — e+7®° and other decays [3] allowed instead. The 
decay of the down quark in the proton, generating wud > uuu+e7 + ib, 
nevertheless, may be considered. Suppose that the distance between the each 
pair of quarks is dgg. Since 


4 e? o+ ie 2 & = 
the location of equidistant quarks in the proton does not effect the binding 
energy. However, now let the d quark decay to a u quark, electron and an 
anti-eletrcon neutrino at the fixed location in the baryon. A reaction d > 
u+W- > ute +, together with p > (uuu) +e7 + may be considered 
with M(uuu)wm, initially. By energy and momentum conservation 
(7) Eq = Ey + Ew- 

Da = Pu + Pw- 


(6) Eetec(uud) = 
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A d quark can decay into an u quark and a W~ boson with a non-zero scat- 
tering angle. The momentum transfer may be sufficient to deflect the entire 
uuu state simultaneously in the same direction. Otherwise, this composite 
state would disintegrate with the third u quark moving away from the other 
two quarks. In the latter circumstatnce, the alignment of the momentum vec- 
tors would require either the W~ boson to move in the same or the opposite 
direction as the u quark. The conservation laws then become either 


(8) yamac? = tuTMyc? + YWw-mMy- c 
YdMdUd = YuMuYu + Ww-My-vy- = 0 

and 
(9) yw- (va Vw-) _ Mu 

Yu(Vu — va) my- 
when the W~ boson moves in the same direction as the wu quark or 
(10) yamae? = Yume? + yw-myp-e 

YdMNdVd = YuMulu — Yw-My- vpy- 

and 
(1 1) Yw- (va ate Uw ? Ma | 


YulVu — vd) = My- 
if the W~ boson moves in the opposite direction to the u quarks. Under both 
conditions, there exists a range of velocities vg, vu > vyw- and vg < vy < 


m2\ 2, m9 2 een . ; 
(1 - mt) ce + FHV; which include solutions to the above relations. 
d 


Therefore, even under the constraint of the alignment of momentum vectors 
in the reaction p > (uuu) +W- -+ (uuu) +e” + 0°, the kinematics allow the 
transformation to the new composite state. The influx of a positive electostatic 
potential energy generates an increase in the mass of the (uuu) state. 
2.0395145184 x 10-°eVm 

daa : 
Given that ma++ = 1272 MeV/c’, equality of the shift in the electrostatic 
potential energy with the difference in the rest mass energies 
(13) Mat+e2 — Myc? = 333.5MeV 


requires 


(12) 


e Sele e? _ 2.0395145184 x 10~°eVm 


4 
(14) Fetec(uwu) = 5 


and 
(15) dqq = 6.115485 x 107'8m. 


- 301 - 


If the only force at this distance is electostatic in origin, ther would be nothing 
preventing the u quarks moving radially outwards and the composite state dis- 
integrating. Consequently, the onset of the strong interactions is necessary to 
restore the constitution of the eventual (uuu) state, which then as a mass in- 
creased by 333.5‘MeV/c?. Equilibrium would be reached at a distance greater 
than or equal to 6.115485 x 107!8m. 

The kinematical conditions also may be given when the byproduct has a 
greater rest mass. The A** baryon is both heavier than the proton and has 
angular momentum 3, Given the decay d > u+W7~ + u+e7 + i, the 
reaction p -> Att + W7 may be considered at non-zero momentum since 
charge, isospin and angular momentum are conserved. Again, the momentum 
vector of Att wil be selected to be aligned along the direction of the motion 
of the proton, and similarly, py- || @. By the conservation of energy and 
momentum, 


(16) Ey = Batt + Ey - 
Dp = Pat++ + Pw-- 
Suppose that p,++ and py- are pointing in the same direction. Then 
(17) [Pol = a++| + Pw-| 
From the two equations 


(18) Wipe? = YareMar+e? +7 -my-c 


Vp pp = YAt+Mat+Vat+ + IYw-My-Uy-> 


(19) Yat+™Ma++ (Up — vat+) + Yw-My- (vp — vw-) = 0. 
Since mp < Mg++, My-, Up > Vat+, Uw-- The above relation produces a 
contradiction because vp — vA+4 > 0 and vpuy- > 0. 


Now suppose that p~a++ and py- point in opposite directions. It follows 
that 


(20) [P| = Wa++| — Pw-l- 
and ; 
(21) Yp™pUp = YAt+Mat+Vat+t — Yw-My-vypy-: 


Combining this equation with the energy conservation relation gives 
(22) VAt+MAt+ (Up — Vatt+) + YWw-My-(vp + vw-) = 0. 


A contradiction results if vp — va++ > 0. Let vat++ > vp. Then ya++ > Yp 
and 


(23) YAt+Mat++ + Yw-My- > YpMp, 
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which contradicts the conservation of energy. Therefore, the reaction p + 
Att + W- cannot occur when the momentum vectors are aligned. These 
contradictions can be circumvented by decreasing the components of the mo- 
mentum along the direction of the d quark motion and allowing the projection 
along of the u quark and W boson momenta along the perpendicular to cancel. 


3. Scaling of the Masses of the Leptons 


It is evident from beta decay that the down quark must contain an electron 
and the radius would be much larger than that of the up quark and approx- 
imately equal to 3 x 107!8m. The difference in the masses of the up and 
down quarks may be attributed to nonperturbative effects [4]. The electron 
and the anti-neutrino then can be traced to characteristics of infinite-genus 
surfaces. Since the anti-neutrino is customarily regarded as a point particle, 
the identification with a point source of non-zero harmonic measure on the 
ideal boundary is consistent. 

The relative size of the electron and neutrino may be described. Point 
sources on the ideal boundary represent ends of the Riemann surface, and with 
a boundary of the end being diffeomorphic to a circle, the electron and neutrino 
shall be viewed as circular strings. A spherical model of these particles would 
be generated from a rapid revolution of the string about its axis. With a 
uniform linear density, the greater extent yields the larger mass. By contrast 
with the Lorentz formula, the mass increases with the size of the circular arc. 

The proportionality of the ratios of radii with the ratios of the masses yields 


(24) Tue = ae Te. 
Me 


The electron neutrino mass is measured to be 0.007614 eV and 


(25) we — 6.71132 x 10”. 


Ve 

It follows that 

eS ae 3 x 1071°m_ 
& > 6. 7lia2 9¢ 10% 

There is a scaling of the masses of the leptons from the first to second gener- 

ations and second to third generations. Let (Be =a and (= = b. By the 

Koide relation [7], 


(26) = 4.47 x 10774. 


(27) Lath =F ta+b), 
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A rearrangement of terms gives the quadratic equation 


(28) b? —4(1+.a)b+ (1—4a+a7) =0. 

Of the two solutions to this equation, 

(29) 2(1 + a) + V3(a + 2), 

the lesser root is excluded because m; > my. The larger root is 
(30) b= (2+ V3)a+2(1 + V3). 


The values a * 14.5, b = 59.578 generate the experimental charged lepton 
mass ratios. 

A theoretical basis for the mass scales of the leptons in the fermion gener- 
ations requires an explanation of the relation between the supercharges and 
the number 14.5. The example of a massive representation of N = 2 Poincare 
superalgebra 


(31) {Qo Qab} = 20" aa Puds 
{Q2, Q3} = 2eapeZ 
3 
{Qh Qh = 22€ 4 5€adZ 


furnishes the anticommutators 


(32) {aq, 05} = bap(M + V2Z) 
{ba, bb} = bap(M — V22), 

where 

(33) aa = 5 (Qh + <op(Qh)') 


(= : (24 - éap(Q3)") 


The central charge then can decrease the trace of the anticommutators of 
supercharges. For example, when the number of spinor indices is increased 
from 2 to 4 in an N = 4 superalgebra, tracing of a set of anticommutators 
can give a factor of 16(M — 3/2Z). Setting the parameter Z of internal 
rotations commuting with the supersymmetry generators equal to KM , the 
factor becomes (16 — /2)M. 

The ratio [Be can be set equal to 16 — 2 only if the supercharges rather 


than the masses are scaled by this factor. Therefore, a Za graded symmetry 
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that squares to a supersymmetry is necessary. One of the anticommutation 
relations of the Z4 graded algebra would have the form 


(34) {Q1°a7, Q% tea t= (1)ay(Qu)ga» 


where two indices (ad) correspond to a single four-component, index. Since 
&Y can cover only four indices, the elements of (y")ay belong to the diagonal. 
Then (1")a4(Qu)s gj can be written as (1) eo (Qy) 36°", where repeated indices 
are summed and Qp is set equal to the supercharge @. The scales of the 
supercharges may be arranged to increase by 16 — V2 from the first to the 
second generation, while there would be an overall scaling of (16 — /2)? for 
the masses. 

The ratio of the lifetimes of the muon and 7 lepton is known to be a constant 
multiple of the ratio of the fifth powers of the masses. By contrast, the Lorentz 
formula requires the ratio of the radii to be proportional to the inverse of the 
ratio of the masses. It follows that 


(35 ry = were = 1.4509 x 10718m 
7) Mu 


rp = ere, = 8.628 x 10-m 
Mr 


Corrections resulting from the weak interactions may be deduced from the 


width of the resonances [ = oe Setting Hp equal to Mpc = 105.658 MeV 
or m,c? = 1777 MeV, the shifts in the mass are found to be 


(36) Am, = 10-“MeV/¢? 
Am, = 1.61 x 107°MeV/c* 


and the effect on the masses and radii is negligible. 

There have been various figures given for the masses of the neutrinos for each 
flavour and the mass eigenstates [8]. A comparison with the Koide relation is 
compatible with a normal hierarchy. The neutrino mass matrix [1] requires a 
mixing angle [2] larger than the Cabibbo angle. Generally, the masses derived 
from the mass matrix and the eigenstates yield an inverted hierarchy with 
values of my, and m,, which are generally not sufficiently large [9] 


4, Quarks in the Second and Third Fermion Generations 


The mass of the strange quark exceeds that of two @ anti-quarks in chiral 
perturbation theory by less than 1 MeV/ c?, when the value mae = 243.15 M eV/ C2 
is used. Combining two @ anti-quarks with a positron produces a composite 
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system of charge —fe with mass 
(37) IMac + Me+ = 486.811 MeV. 


Beginning with a dd state, a charge of Re would have to be added without 
introducing a d anti-quark. Instead, contact with a d anti-quark must result in 
the transfer of the charge ze and a recoil of the anti-quark for this configuration 
to model a particle of charge —fe and mass nearly equal to msg. 

Both composite systems, however, do not satisfy a mandatory condition 
on the representation theory of SU(3),. Quarks transform under the 3 rep- 
resentation and anti-quarks belong to the 3 representation of SU(3). Since 
the positron is an SU(3) singlet state, it is not tenable for tet to transform 
under this group identically to the strange quark s. Similarly, the dd(3) state 
transforms under the 3 x 3 = 8 @ 1 representation, which is also different. 

An alternative configuration consistent with representation theory of SU(3) 
would be the a we~ state together with four-infinite genus components with a 
mass 
(38) 

Mue + Me- + Amuco = 243.15 MeV/c? + 0.511 MeV/c? + 4(61.85 MeV/c?) 
= 491.061 MeV/c’. 


The sum of the radii of the chiral up quark and the boundary of the infinite- 
genus components equals 


1, 
(39) Tuc t 4fyco = Tue + sare = 1,45244618676 x 107m 
which is much too large for the strange quark to fit inside a baryon from the 
decuplet that includes the Q~. Instead, the boundary of the infinite-genus 
surface must be wrapped around the fundamental constituent up quark. Its 
final radius is measured to be 5.95268 x 107!9m. Setting 


N 
(40) 2m 5° (2.808563687 x 107m + nAx) = 2n(1.45244618676 x 10-'3m), 
n=0 
it follows that 
(41) N = 3.31562 x 10° 


Aa = 9.4827406 x 107m. 
The electrostatic potential energy between the u quark and the electron 
De 36? __ 1.0705775916 x 10-%eVm 
~B4meodue- ue 
must be included in the computation. When this electroatic energy is dis- 
charged into the bound state during the process of fusion at a distance of 


aa) ve 
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due~ = 3x1071%m, VS = 3.598591972 MeV yielding ms = 491.061 MeV/c?— 
3.598591972 MeV/c? = 487.462408028 MeV/c? which belongs to the range of 
experimental error. 

The charmed quark may be described as a complex consisting of dd, u and 
dd, with one u, two d and two d infinite-genus components. Then the mass 
would be 


(43) Me = Mue + WMde + 2M, + Muoo + WwNdoo + 2AMNgoo 
= 5(243.15 MeV/c”) + 4(64.85 MeV/c”) + 61.85 MeV/c” 
= 1543 MeV/c? 


which is consistent with the mass of the J/% particle. The electrostatic 


potential energy of the complex describing the c quarks is Vj, = ye + 
ove. + Vee. Vb OVI a Vow ‘, When the u quark is placed cen- 


trally in a square with d and d quarks at the corners, V%, + ue = 0 and 


i elec ele 
Vad + Ve + 2Vaanor, = 0. The remainder is 


dddiag. 2 e? 3.598591972 x 107eVm 
ae) 2Velec | = 'QAveiaae 8 ae = eS 
TEOGdddiag. dddiag. 


When the distance between the d and d quarks across the diagonal of the 
square array is 8.3688185395 x 10-!8m, the mass of the charmed quark is 
reduced to 1.5 GeV/c?. 

The quarks in the third fermion generation require many more components. 
The YT resonance at an energy of 9.46 GeV is a17 bb state if my, = 4.73 GeV /c?. 
The following configuration 


U U U U U WU 


U a u uu a 
has a mass equal to 


(45) 6mue + 6Mage + 2Mde + MG + 6Muco + BMaco + 2Mdoo + MGog 
= 4584 MeV/c? 


The difference from the experimental value of the 6 quark mass can be ex- 
plained through the addition of infinite-genus components neighbouring quarks 
and anti-quarks and a reduction as a result of the electrostatic potential en- 
ergy. If the d and d quarks are placed symmetrically between the uu and 
uu pairs, the electrostatic energies involving the u and & quarks will cancel 
with the exception of neighbouring the diagonal interaction between u and u 
quarks. When the planar array is replaced by a connected three-dimensional 
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configuration, there is a negative electrostatic potential energy between the 
uu in the third complex and the uu pair in the first complex. If the three 
complex are arranged in an equilateral triangular form, the distances between 
the central d and d quarks will become equal. 

Then the electrostatic potential energy of the entire supercomplex would be 


(46) 12Vua + Vaa + Vag + 2Vaq + Vaa + Vag = 12Vua + 2Viy 


These electrostatic energies equal 


(47) puadiag. = # eF _—7-197183944 x 107MeVm 
eke 9 Ato dyadiag. Ayadiag. 
yaa 1 6? _ _1.799295086 x 10-eVm 
elec 94neodaz dag 


The length of the diagonal from the u to & quark is V2 times that of the length 
of side of a single complex. The distance between the d and d quarks would 
equal the length of this side. 

an the distance between the surfaces of the d and d quarks is 6.335412 x 
10-*"m, 


(48) ui — —80.3290105 MeV 
Vad, = —28.4 MeV 


The net electrostatic potential energy is 
(49) 12Vua + 2V43 = —1020.748126 MeV. 


The strong interaction potential would contribute significantly to the energy 
unless it drops rapidly to zero beyond a range between 107m and 107!7m. 
Since 6mMucot6Macot2MdootMeog = 12(61.85 MeV/c?) +3(64.85 MeV/c?) = 
936.75 MeV/c”, the mass of the entire supercomplex of 15 elements is 


(50) 4584 MeV/c? + 936.75 MeV/c? — 1020.748126 MeV/c? 
= 4500.001874 MeV/c’. 


It may be noted that the sum of the masses of the supercomplexes and 
the boundary components gives 4584 MeV/c? + 936.75 MeV/c?. Therefore, 
the YT resonance would be created by a fusion of the b and 6 quarks with 
the electrostatic potential energy of a single complex approximately equal to 
—184.8166 MeV, which can be achieved by increasing the distances between 
the d and d quarks to 8.1781346 x 107!8m. 
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The top quark is much heavier with a mass of 175 GeV/c”. Let n be number 
of complexes of the form 


d d d d 


Each complex has a chiral mass of 5(243.15 MeV/c?) = 1215.75 MeV/c?. 
Given that there are infinite-genus components of surfaces emanating from 
each of the chiral quarks, the entire mass contributed by n complexes 

(51) 

n(1215.75 MeV/c*)+4n(64.85 MeV/c”) +n(61.85 MeV/c?) = n(1537 MeV/c”). 


Setting n equal to 113, it would equal 173.681 GeV/c”. Given that there are 
two other infinite-genus ends originating from the central 113 w and @ quarks, 
their connections generate another 


(52) 226(61.85 MeV/c?) = 13.9781 GeV/c’. 


Adding the two masses gives an estimate in excess of m; by 12.6591 GeV/c’. 
The remainder of the supercomplex would close to form an electrically neutral 
structure. 

There is a reduction in energy resulting from the electrostatic forces between 
the a three-dimensional polygonal configuration with alternating central u and 
@ quarks. With 57 wu quarks and 56 &@ quarks in the complex, the electrostatic 
potential energies amongst the central quarks and anti-quarks cancel except 
for 56 x 2 V%#. The diagonal interactions yield 113 x 4 Vie 


elec’ 


Setting dyg = 1.09104255722 x 1071"m and dag = 1.54296713803 x 10-!"m, 
(53) wie, = —65.9661146814 MeV 
V4, = —11.6612722439 MeV 


elec 


The mass is reduced by 12.6590998986 GeV to m: = 175 GeV/c”. 


5. Conclusion 


The mass spectrum of hadrons has been successfully described through the 
quark model and the spin coupling formula. Further regularities in the inte- 
ger multiple rule, Regge trajectories and the Koide relation have resulted in 
substantial understanding of the elementary particle interactions. The pat- 
terns in the fermion multiplets of the standard model have yet to be entirely 
elucidated. The quarks are found to have a mass generated by chiral pertur- 
bation theory and a nonperturbative portion that may be attributed to the 
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infinite-genus term in a string theoretic version of a series expansion of expec- 
tation values of products of field strengths. Many of the conventional methods 
may be adapted to the chiral mass of the up quark, which is 243.15 MeV/c’. 
This value is sufficient to establish an inverse proportionality of the ratio of the 
masses of the up quark and the electron to the ratio of the radii by the Lorentz 
formula. A similar calculation yields the radii of the muon and 7 lepton. Since 
the beta decay requires the reaction d — u+e~ + %, and the difference be- 
tween the masses of the up and down quarks is a nonperturbative effect, the 
electron and the anti~neutrino are identified with characteristics of the end of 
an infinite-genus surface. By constrast with the Lorentz formula, the extent 
of the boundary circle representing the two particles is proportional to the 
mass, given a uniform linear density. Then the radius of the anti-neutrino is 
computed to be eight orders of magnitude less than that of the electron. The 
analogue of the beta decay for the proton, converting a uud to a uuu state, is 
found to be kinematically impossible when its momentum is zero or aligned 
with the momentum vectors of the byproducts in the intermeidate channel as 
a result of the larger rest mass of the A** baryon. Therefore, it is verified that 
a proton at rest is stable within the standard model, and a similar conclusion 
is valid when it has momentum unless a composite wuu state of lesser rest 
mass is created or the momentum transfer allows a non-zero scattering angle. 
The Koide relation for the masses of the charged leptons has a theoretical 
explanation in terms of an anticommutation relation for the supercharges and - 
an approximation for disparate masses. Definining two variables a and b equal 


to are and ae respectively, a quadratic equation for bin terms of a may be 


derived, with the solution =(2+/3)a+2(1+[V3). The derivation of the scaling 
of the first generation, a, requires the introduction of a coefficient in the trace 
of the anticommutators of the form 16 — V2, which is possible with central 
charges in’ the algebra. A consideration of four-dimensional superalgebras 
reveals that this factor can be introduced for masses. However, consistency 
with the scaling a requires that it is present for supercharges, which requires a 
Za graded supersymmetry. The specific form of each of the anticommutation 
relations, and the central charge, would be a necessary for this theoretical 
explanation to be valid. 

Given the infinite-genus components representing the nonperturbative con- 
tributions to the u and d masses, various models may be constructed for the 
quarks in the second and third fermion generations. Group theory of the 
SU(3)- makes it mandatory that the configurations transform under the 3 
representation. An equivalence between the strange quark and a combination 
of the up quark and the electron can be completed when four-infinite-genus 
components are added, the electron is scattered and the ideal boundary is 
wrapped around the chiral quark. Complexes of up and down quarks and the 
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anti-quarks with five elements yield a composite system that satisfies this con- 
dition, with the correct net charge. Then a sequence of identical copies of the 
complex and the antiparticle version would comprise a much larger structure 
which models the more massive quarks. The central quarks are attached to 
two infinite-genus components in the third generation, and there is a reduction 
resulting from the electrostatic potential energy. By this technique, the ec, b 
and ¢ quarks can be described by 1, 3 and 113 complexes. This sequence must 
be explained for a rigorous theory of the structures of the quarks. 
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Abstract 

It has been shown recently that superstring states are organized into supermultiplets 
for the massless fundamental modes and solitons representing excited states. Summa- 
tion of superstring amplitudes yields an expression that has poles when the momenta 
of the external states have specific values. It is established that only the residues and 
not the mass can be shifted by higher genus terms in the superstring perturbation se- 
ries. The masses of elementary particles will not be generated by radiative corrections 
in superstring theory. Quantization of the energy levels does occur for the excited 
modes and the solitons. The superstring must be replaced by a hadronic string to 
achieve a phenomenologically viable theory of the particle spectrum. A. linear spacing 
for the masses can be derived from the pion bubble model. The vibrational modes of 
the string between the quark and anti-quark in the pion then generates square-root 
variations from the integer multiple rule for the excited states of hadrons. 
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1. Introduction 


While solitons described by extreme black hole geometries might be used to 
representing elementary particles, excited states generally have a Planck mass, 
and it is useful to restrict consideration to the massless sector. The feasibility 
of such a truncation is raised by connections between Type IJA string theory 
and eleven-dimensional supergravity [23] [24]. 

An alternative approach of the derivation of the elementary particle spec- 
trum would follow from characteristics of superstring amplitudes with mass- 
less external states. It may be recalled that the energy levels of string states 
are defined by a Fourier mode expansion and expectation values of of oper- 
ators satisfying a Virasoro algebra. A quantization of energy is introduced. 
The superstring amplitudes will have poles at the on-shell masses of excited 
modes defined by the external vertex operators. It is analogous to the poles 
of the Koba-Nielsen amplitudes in bosonic string scattering. Corrections to 
the masses at sufficiently high order are prevented, however, by nonrenormal- 
ization theorems. Given the quantization of on-shell energies, it would follow 
that there is a theoretical basis for the integer-multiple rule for the elementary 
particle spectrum in a superstring model. 

An evaluation of resonance poles in superstring correlation functions yields 
a series of corrections to the residues without affecting the masses. The ab- 
sence of a mass renormalization would be expected from the finiteness of the 
superstring perturbation series at each order and summed over the genus for 
a range of the string coupling. Therefore, higher-genus terms cannot shift the 
massless states to the masses of elementary particles. 

The predicted values of the masses based on the excited states of the super- 
string would not be considered to be phenomenologically viable. It is necessary 
to replace the superstring by a hadronic string. If the hadronic string model 
has supersymmetry, this symmetry would be present in the spectrum [10], 
which may be valid for two multiplets with a tetraquark resonance. By con- 
trast, a nonsupersymmetric hadronic string model can generate the known 
spectrum, but the quantum theory has problems of consistency similar to that 
of the bosonic string. 

The integer multiple rule, nevertheless, characterizes a branch of both the 
meson and baryon resonances. An explanation through pion bubble diagrams 
[11] is compatible with validity of the rule for the 7 branch. This method 
then may be combined with the hadronic string model to derive a realistic 
spectrum and a square root variation from the integer multiple rule for the 
excited states. 
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2. Scattering Amplitudes 


The poles of string amplitudes have integer spacings. The scattering am- 


plitude with N vertex operators would be 


(1) Doe ee? dure f, llé ay V9(%i) (Vi (21). Vur(Zn), ). 


Xg t=1 


where V;(z%) is the ith vertex operator, duwp is the Weil-Petersson measure 
on the moduli space M, of complex metrics on a Riemann surface of genus g 
and Kg: is the string coupling. The path integral of the time ordered product 
would be evaluated through the Wick contractions of the operators, which 
yields 


(2) 


( 


0 (1 + tp} -X(z1) = wr ‘ X(z))? _ a (Pr: X(z))? + -) 
(1 + tpg + X(z2) ~ si? - X(z2))? - ap." - X(z2))* +. .) 


(1 + ipy + X(zn) — aN -X(zn))? — a(n -X(zn))? + -) 0) 


= (0|0) ie (Olp: » X(z)|0) — $5 Olpitpy Xu (22) Xv (zs) |0) 
t<j 
+ DY (lpttn Xu (2) Xv(24)10) Olek? Xp (zx) Xo()10) 
eae 


on 20 pt'py pi pF (0|X u(zé)X p(zj)|0)(0| Xv (a) Xo(zs)I0) 
1<j 
+ D0 mi prope (0|X (21) Xv(24)|0)(0|Xp(zj) Xo (zx) 0) 
ea 


“oF mana 8 > _ Pipi nh § 975 (0|X (zi) Xp(zs)10) (O| Xv (24) Xo(z4)10) (0X (%)X-(z)|0) 
i<j 


Cr 2” Pipi pp pSpepE (O1X p (zi) Xa(zj)10) (OLX (zi) Xo(z4)]0) 
i<j<k 


(0|Xp(z5)Xr(ze)I0) 
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(3) 


— Gin? Do whens | (01x (2) Xe es) O10Xa (2s) Xo (4)10) 
; i<j<k<e 


(0|X p(z%) Xr(ze)|0) 
wags aseb ap tyteea sacccanesesdenese + (0|X.(24)Xp(2x)|0) (0| XL (23) Xo(Ze)|0) (O|Xy (23). X+(ze)|0) 
+ (0|Xp(zi)Xr (ze) |0) (O| Xv (25) X p( 2%) 10) (OLX (23) Xo (ze) |0) 
.. + correlation functions with coincident fields 


+ expectation values of normal ordered products of operators 
+... 


The general combination is 


(4) (Diy * pH)" «(Din Pim) (G (eins ZA)" ++ (G (Zinn Zin) ™ 
Tit. tlm =n. 


Suppose that 71, 71, 22, ja, «> im; jm are different. The combinatorial factor 
from the expansion of the correlation function is 


(5) 


From the expansion of the exponential exp (- icy PE PY (OX (zi) X (z;)10)) 
is 


1 n 1 
(6) ni (,, adi vn) ~ ryhtml 


4 


If some of the indices 71, 71, 72, jo, .-.-4m jm coincide, an alternate labelling 

k1, ..., ke of the expression is 

(7) | 

(Pk1‘Pko)”* (Pko"Pka)”?+++(Pheig Phiig tr) "0 ++(Phe1 Phe)” G(Zhy » Zhe)? --G(Zhep_15 Zhe)” 


where pp, appears Thy = 1ky,3;, + + r1,81,t, = Thi times, ..., Dp, appears 
The = Vheso1 + ++ + The,set) = The times. The combinatorial factor from the 
expansion of the correlation function is 
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See hs i‘ he pt 
PhyteThe! \ VkyystaeT ki, 81,4, Thesde,.°7*T he,S0,t, 


1 
= eo 
ky ,814 They 81,04 tT hee ,8¢1 tT ke 82,t9 H 


1 
r!.-Tm! ; 


The expansion of the exponential again gives ok Therefore, the coeffi- 
cients are equal, and expression for the correlation function is given by 


(9) ] [ewe (pi: ri (a, 25) 

i<j 
with G(2z;, z3)n"” = (O|X"(z), X”(z;)|0). and not including the normalization 
of the Green function. The Green function for the free scalar field on the 
complex plan would satisfy 
(10) d0G(x%, 23) = 5(% — z;) 
The solution is 


1 
(11) G(%, zj) = ain \24 — | 


since 


(12) 08 5(In 2+ In 2) = 3 (=) +2 (=) 


and 
(13) -5 [ eex[0(2) +9(2)]=-42 


On a Riemann surface, two delta function sources are necessary for consistency 
with Stokes’ theorem, and a summation over image charges under the action 
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of a uniformizing Schottky group gives 


(14) 


au zp — VozR 2Q ~ Vazs 
Gas(P, R) - Tv le 2Q Vazg 2Q VoZR 


Qa 


— EY Re fomer) —am(za)} Em rYahRefen(en) ~ ot25)} 


myn=l1 


Un(z) = d ln (Zee a Un(Z) ~ Un(Imz) + 2ritmn, 


where R and S represent the sources, 7 is the period matrix, and {vp, n = 
1,,,-g} is the set of integrals of canonical basis of homolorphic differentials and 
the sum over elements {V,} over the free group Coe separ ye does not 


include those products of generators that have T+! as the right-most factor. 
The final expression 


(15) II em ZPiPiG (224) 
i<j 


then would be required to obtain an integral of the form 


(16) 


— P14 
t(28, 29, 28) | Exaile — 20) TT 2 
A o£] 


* Tex Le PA Reef Um(z4) — Um(22))(Lm 1)~1 Re(vn(z,) — Un(28))} 


—PLP4 
— Hest = es 2 


x (similar factors with 22 - 23, 23) Pi —> Po, 3). 
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The four-point finite-genus amplitude, for example, then may be expressed 
in the Schottky group coordinate parameterization of moduli space by 


(17) 
[tele AE fe BE — BPP Oe, 2) 
20 . ao « _ \l(s _ 50\m 
¥ yh de Pe a (2 — 21)"(2 — 21) 
DT aca 8 42 Ale PE YY SE 
— Pap. 7 
HI {|z — 29|- WP |z — BIE O(z, 2} neo 
on (z — 29)*(z — zo)™ 
+f I dé diz —z 27 2 O11— i tas (2 = 2a)"(% — 3)” 
F lea <A | 3| | 3| =D film! 


POM — AOE Ne AE Be hae 


2m PPA (z — 29)*(z — z)™ 
+f i dO dlz — 22| |z— 28}t- (eiita) (2 2a)” 
x AISA | 3| | 3| > E: lm! 


Phat 
2 ea 7 
Hz — O|- a | 2 — Bl 2 B(z, Dhemed 
+ finite 
—OP4 +on+2 or 
= apa ea) {lz — 28) 2 — 28) Bz, 2) } 0 


An Pt +2n+2 


_PLPa _Pa-p e 
+ Lam FTO are 5 Tape OO - Bye 2 — Bl 8 (2, Z)} 20 


—PaP4 4 on+2 pupa 
# > Sear” {lz — aI 2 — QI Oz, 2) baa 


+ finite 


where the domains of the integrals are divided into {z||z—z?| < A}UF,/{z||z— 
2?| < A}, i = 1,2,3 in the Schottky covering surface CUoo/U?_, (Dr,,UD7-1), 
where ®(z,Z) is a regular function. The poles in the amplitude occur at 
s, t, w = 8(n— 1), beginning with the tachyon and progressing to excited 
levels. When n = 0, the poles occur at s, ¢, uw = —8, which would describe an 
intermediate tachyon resonance. Another linear dependence, such as 4n — 12, 
which would arise if the exponents included —p; - pg without the factor of 
i, yields intermediate states with squared center-of-mass energies equal to - 
12 and, after iteration of the scattering calculation, a set of resonances with 
arbitrarily large negative energies. Consequently, the factor of $ would be 
present in the factor (15), and the tachyon vertex operator actually must have 
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the form eiV GPX). The spacings are determined by the tension of the string. 
Therefore, the energies can range from Planck to hadronic scales. 
An N-point correlation function in ah theory would be equal to 


(18) (Yo. Vy) = oe? & duwp |] de(Va(er) « Vielen) 


g=0 t==1 


Expanding about a particle resonance 


(19) ee ee 


g=0 

It may be recalled that the upper bound for the amplitude at genus g is 
given by the maximum over the interior and the compactification divisor. 
The integral over Schottky group parameters over the interior, which is the 
complement of the neighbourhood of the divisor, is less than z, where the 
quotient by g! results from the restriction to the fundamental domain of the 
supermodular group [6]. Summing over the interior and the components of 
the compactification divisor describing two components, the bound for the 
sum over the amplitudes would be 


mos 5 Bema Sst] 
ge 91 +92=9 n! ge a: gitog=o  91' 92" 
191 ,‘9g<n-1 1$91; 99<n-1 


Let 


11 
(21) X= YD == 
g1t+92=9 gu} ga! 

1<91,92S59~-1 


Since 


2 


(22) etite2 yy Gm (wi + %2)? + ae = 3) » aft af? 


et? 
g=0 g=0 91+92=9 gu! 92: 
[oo] co 
(23) e=yF =~ > ala = S2y 
g=0 g! g=0 gitga=g"t' 77" g=0 
and 
29 
(24) Pag = oy 


g! 


Then 


1 1 1 1 29 —2 . 
(25) 2. a Oa gio g21. 
aiteg=g 91" 92" gS F g 


1$91,:92Sn-1 


The magnitude of the amplitude is less than 
gad. _ 
(26) K29-? | Bis Ei + ee = K°9-? |Bl9 ae 
g! g! g! 


Then the summation of the factorizations over the particle resonance is bounded 
by 


(27) ee ae 


—m 
g=0 


The series converges for all values of the coupling; and the mass of the reso- 
nance remains invariant. 

Infrared divergences resulting from propagstion of tachyon states in bosonic 
string are cancelled by supersymmetry. Then there are massless states propa- 
gating along thin tubes created by pinching the separating cycle between the 
genus k and genus g = k components of a surface of genus g. After these 
resonances are excited states of the string. The problem with this model is 
that the excited states of a Planck size string have Planck masses, which is 
too large for the elementary particle spectrum. The viability of a model of 
the particle spectrum based on superstring amplitudes would depend on ex- 
istence of corrections to the mass of the resonance. It is evident from the 
above series that there is no mass renormalization in superstring perturbation 
theory. Only the residues, representing couplings to different channels in the 
scattering process [3], are altered by radiative corrections in the superstring 
perturbation series. 

The nonperturbative contribution to the amplitude arises from the infinite- 
genus boundary components. The amplitude given by the integral over the 
set of effectively closed surfaces with zero harmonic measure would have the 
same form as finite-genus amplitudes and not affect the masses of resonances. 
However, there exists a class of surfaces beyond this category of Og surfaces 
that represent the nonperturbative phenomena in the theory. A nonvanish- 
ing capacity of the ideal boundary yield a non-zero potential for the surface 
through the integration of the normal derivative of the harmonic measure 
[16]. The potential energy can be set equal to the energy of the infinite-genus 
boundary component, When there is no other dynamics, this energy would 
be proportional to a non-zero rest mass. 
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It is known also that the masses of leptons may be estimated through ra- 
diative corrections [20]. If the Planck strings are replaced by hadronic string 


2 
with tensions that are multiplied by a factor (n=) , the masses of the ex- 


cited states become more realistic. There is still a large difference between the 
masses of the hadrons and the leptons. 

The formula for the mass of the electron in a quantum theory of radiating 
oscillators is precise, however, and it is indicative of the necessity of methods 
which are valid for the finite size. The equation for a radiating electron in an 
electromagnetic field is 


, 2 &a 
(28) Max — 303 de =e 
with E = E,,cos(wt) has the solution 
(29) x(t) = — e005 Cycos(wt — G,) 
where tan ¢, = =o The Hamiltonian may be given by 


(30) A! = 4c? Sw,” b» eiBisin SsicosCgsin(T'si + C54)” 
38 t 
+ oy excos C5408 (I's: + C5s))” 
i 


and #,; is the angle between the wave with index s and the velocity Bi, while 
C4 = are tan (-#), with wo = Smee and wi; = Ws(1 — Bicos B.i)(1 — 62)-3, 


and ['; = coericus tat 4. [17]. For a sels particle at rest, 8; = 0 and C/, = ¢s. 
It follows that 


(31) E} = dnc? e? Dawe cos*¢,cos* (T's; + Cs). 
Replacing the sums by an integral, 
(32) Fo = dmc? i / dw (1422) cost AC) 
i w toned we Brea, 
and setting (cos?(I's; + ¢s)) = $s 
2 poo 2\ 1 2 2 3 m.c2 

0) iil = cio | Smic _ 3 

(33) (Ez) = = | (1 + =) dw 7 mic’. 


For the interaction between two charges, it may be noted that cos(I'4; + 
¢s)cos(P's;+¢s) = $cos(Pss + T's; +25) +§ cos(I's;—I's;) and, when cosy, = 0, 
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1 _ 1 Ws(%3—1T;)cos 9. _ 1 Wslij _ = 
3cos(l'sirsj) = 5c08 Geacanee s—"i) +) = jcos (74). Let ny =r, & = uao 


and xo = =, where rp = eer. Then 
(34) 
e2 2 ee Cc 2 1 
E}.= ef du-aysin(uao)(1 + u?)~* = <2 [* dusin(uxp)(1 + u?)~*. 
CT 0 rT 


The evaluation of the sum of oscillator contributions is the equivalent of a 
one-loop calculation of the self-energy of the electron. There also would be 
higher-order contributions in the diagrammatic expansion of quantum field 
theory. If the action may be included in a supersymmetric model, a nonrenor- 
malization theorem may be used to render corrections to be vanishing. When 
supersymmetry is not present, the magnitude of any higher-order terms would 
have to be bounded. 

A classical theory with form factors would begin with, the replacement of 
the current by e [°° v,F(R*) in the source term f j“(z)Ay(zx)d*z, where 


s K4 
F(R) = gy fe" 8 9(K*)d' K, g(K’) = mReTK’ and terms in the elec- 
tron equation include (a? +e ? Sa.) by, — 3e2(%, + 07v,) + ..., where Sa = 
foe, H(o?)do and H(R?) = —4n f PU) ei ‘Rat Kk [4][18].The electromag- 
netic contribution to the energy is equal to zen if 


(35) i he H(o2)d 


ik 


3. Integer Multiple Rule 


The integer multiple rule has been observed amongst leptons and hadrons. 
It may be noted that the formula for the masses of dyons in N = 4 super- 
Yang-Mills theory [ 22] verify the integer multiple rule when either the electric 
or the magnetic charge is significantly larger [5]. The dyons represent nonper- 
turbative states in string theory. 

Nevertheless, shifts in the masses of hadronic states are expected to be 
derived from radiative corrections. By the formula for the bosonic string 
amplitude, it is evident that the poles occur at 


(36) Pitt = mn+2 
or s, £ ,w =8(n—1). Consequently, the squared center-of mass energy scales 


as the square root of n. The path integral formalism is developed to com- 
pute the scattering amplitudes of N massless bosons with vertex operators 
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Vo = Ey (Ox — ip php” )eP® and V_ = u*(p)Sa(z).Sq *(e)etne and 2n mass- 
less fermions with vertex operators V_1 = u® (p)So(z)Sq ?(z)e* with ghost 
charge —5 and Vi = u%(p)Sa(2)S2 e?® with ghost charge }, where Sq and 
Se 2 are spin fields [1]. The fermionic scattering amplitudes 
(37) 
—prpo+Goy Gag ~F-+n+1 


N ror 
Ay = [ dxo(epatin I di;w(z;) S- (w&) 
j=3 


! 
mar ni 


ve [Or (21 ; 29) F(z » 225 23)| lees 


N 
F (2,2, 23) = [] O(a, zj)tar GI F1+F O29, zj)ten GPa Bit 
j=3 


1 21 1 22 zy -1 
len -5 u-5[ wt Dos [ | Kij 


Kij = [J (O(a, 2;)]®O~ By 96-95 
25t<j 


where 7p; and p2 are momenta of gauge bosons, gq, and qa. are fermion charges 
A 
of the spin field, w = wars) , with {w4, A = 1,...,g} is a canonical 


it 
basis of holomorphic one-forms, ge the tangent vector to the geodesic con- 
necting z; and 22 in a normal coordinate expansion on the Riemann surface 
and g; is the ghost charge, has poles at p1 + p2 = Yor * Ya. — i +2+1 f1). 
The squared center of mass energy is scaling linearly with the level number. 
Recalling that, in a center-of-mass frame, s? = (pi +p)", it would equal 


~ 3 
(38) = ph + rh+ 2 doa +2 (m4 9). 
Setting s at the resonances equal to M?, where M is the mass of the inter- 


mediate state, when the vertex operators represent massless states with null 
momentum vectors, 


bat dey 3 
(39) M? = 2, * Fog +2 (n+3). 
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The masses of the charged particles may be deduced from the dyon and 
Lorentz Formulae. It is feasible to demonstrate the integer multiple rule for 
more massive resonances. The use of the nonperturbative formula for the 
electron, however, would require introduce additional charges. 

The intermediate vector boson in the fundamental interactions has integer 
spin to conserve angular momentum in the scattering of fermions with half- 
integer spin. By contrast, Regge amplitudes can have poles at the squared 
mass of baryons with half-integer spin. Baryons are massive compound states 
that have quantum numbers of the resonances on the Regge trajectories. The 
electron is an elementary particle that generally could be interpreted as an 
incoming or outgoing state and not as a resonance produced by the scattering. 

Nevertheless, the formula (39) may be checked in this regard. The values 
of the inner products may be Qo, «Gaz = —3, iy S. Then, when n = 0 and the 
spin field scalar product is 4 

1 3 

(40) Mpno = 2:5 +2-7 =2 
The ratio of the fermion contribution to the entire expression is 3 :2= 
3, which is identical to that of the Lorentz formula for the electromagnetic 
contribution to the mass of the electron. Since the equality (40) describes the 
squared mass instead of the mass, an alternative mechanism for describing a 
spin-5 intermediate state would be necessary. The existence of the poles in the 
scattering amplitude reflects the quadratic form of the propagator for integer 
spin particles. When the intermediate state has half-integer spin, it must be 
modified to be linear in pant Since 
a ee ee 

pe-m ptmp-m 
it follows that if singularity in the scattering amplitude occurs at p = m, the 
pole in the fermion propagator would arise if (p>—m)# = 0, where ~ is the spin- 
3 field, and the product would be Sapom The factor (w£)~P!P2+4ay Jag Ztntl 
in the integrand of the superstring amplitude arises from the representation 
of theta functions producing the expectation value of the product of pairs 
of fermion vertex operators, V_1 and Vi , and the poles in the scattering 


amplitude may be found through an expansion in a Taylor series. The poles 
of the fermion propagator would be expected to arise from a factor with an 
exponent with a linear function of the momentum. The production of an 
intermediate spin-5 particle, however, requires one of the fermionic incoming 
particles to be replaced by a bosonic particle, with the effect of leaving an 
unpaired fermion in the expectation value. The contraction of factors in Vo 
and Vy, 1 would include that of e? and Ox, which yields after integration 
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over the Riemann surface, yields, a linear function of the momentum. It would 
also introduce only a single fermion charge of the spin field qq 

Setting this linear function equal to the mass at the pole of the amplitude 
then would yield the equation, since p? gives 2v, half of the rest of the exponent 
for p would give 


(42) Merin~3 5 lal i (x + i) 


for a fermionic intermediate state. Since |qq| = $; when n = 0, 


win LLB 
(43) M,=0 5 ae ae 
and the contribution of the matter fields would be 3 of the mass of the spin-4 
particle since the remainder would result from the spin field of the superstring. 

The ratio of the masses of the +y branch of the meson and baryon spectrum 
to the pion mass satisfies an integer multiple rule [2][12] [15][19]. Fractional 
powers may result from a variety of potentials in the modeling of hadronic 
spectra [14]. Since the energy levels are defined for the scattering of two string 
states in a fixed channel, it would be best to consider the integer spacing for the 
squares of masses amongst the corrections to the fundamental mass spectrum 
within each particle multiplet. The Gell-Mann-Okubo mass rule is known to 
be linear for fermions and quadratic for mesons [13] [21], which describe bound 
states of two quarks, and therefore, a resonance resulting from the interaction 
of two strings. 

Given a collection of pion bubbles in a meson with a mass larger by a factor 
of uv, the poles in the string amplitudes would yield corrections 5M for each 
pion that are proportional to in. Given k excited energy levels of the open 
string between the quark and anti-quark in each pion, there would be ku 
excited states of the meson. 


4, Conclusion 


The masses of the elementary particles are known to satisfy an integer mul- 
tiple rule. It has been determined that this rule follows from the masses of 
nonperturbative string states which are dyons in N = 4 super- Yang-Mills the- 
ory. The perturbative string amplitudes have poles at values of s, ¢t and u 
that increase linearly with the level number of the excited state of the inter- 
mediate resonance. It follows that this radiative process then yield quadratic 
corrections to a fundamental energy spectrum. The problem of the linear pro- 
portionality of masses with the level number may be resolved for half-integer 


spin particles by transferring the pole of the resonance to the factor —. and 
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integer-spin particles through pion bubbles since the ratio of the mesons in 
the y branch. The validity of the pion bubble model for the y branch follows 
from the decay 7° —} 2y for each pion. The Gell-Mann-Okubo mass rule is 
conjectured to be the primary method for predicting masses differences in the 
hadron multiplets. 

Phenomenologically realistic values of the masses would be derived from 
the hadronic string rather than the superstring. Without supersymmetry, the 
infrared divergences resulting from tachyon states must be removed. There 
exists a method for ensuring finiteness of the string theory having only modular 
invariance and breaking of supersymmetry [9] through a matching of states in 
different sectors. The tachyon then would not generate the problems for such 
a model that occur for the bosonic string. 

Given the meson resonances in the y branch which have masses that are 
integer multiples of m,o or m,-, excited states resulting from the quantized 
vibrational modes of the open string connecting the quark and the anti-quark 
that differ from the mass of the resonance M by 6M, where 6M « i/n. There 
will be kv energy levels when the open string has k excited energy levels. 

The variation from this rule for the v branch is indicative of an extra con- 
tribution to the masses, which may be related to the infinite-genus boundary 
component that is found to be sufficient to explain a nonperturbative contri- 
bution to the 7’ mass [7]. 

Lepton masses would not be predicted by the current string models. A 
recent study of the generation of lepton masses is based on the division algebra 
modules in the spinor space of the standard model [8]. The infinite-genus 
boundary component is not required. However, the interaction between quarks 
and leptons in the weak interactions is indicative of a connection between the 
ideal boundaries of infinite-genus surfaces and the v branch of the meson and 
baryon resonances. The electromagnetic contribution to the mass of a charged 
particle is given by the Lorentz law which depends on its radius. The formula 
for the energies of a set of elementary particle states therefore would include 
the combination of the integer multiple rule and a Lorentz term that represents 
symmetry breaking. 
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Abstract 
The belief that the substratum has a discrete grainy-like 
structure at some small scale follows from the natural scale set by 
of the Planck length and the assumption that the concept of 
quantum correlations preceded any primitive notion of space and 
time. Dealing with the present formalism of quantum theory, it is 
possible to incorporate the notion of a discrete time interval into the 
‘theory’ through :a discrete time difference amended Schrodinger 
equation. One of the consequences of such a theory is that it leads 
to discrete time corrections to spin precession frequencies and spin- 
flip frequencies that depend on the internal structure of composite 
particles. Thus the theory provides us with a potential probe to the 
composite structure of quarks, leptons, gauge bosons and Higgs 
particles. It is also possible that due to the environment a quantum 
system might be subject to Markov environmental influences that 
perturb a quantum system with an assumed Hamiltonian. In this 
» case also it turns out that these Markov effects are sensitive to the 
internal structure of composite particles and thus provide a probe to 
compositeness. The following discussion suggests various models 
that may probe for both discrete time difference effects and Markov 
environmental influences. It is hoped that future precision 
Measurements of spin-precession frequencies, spin flip frequencies as 
well as spectral shifts in atomic systems may lend support to these 
models. 
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1. Introduction 

The present theory of elementary particles has brought with 
it a degree of unification embodied in the presence of a gauge 
symmetry (SU(3),, X SU(2), X U(1),) along with the beautiful idea of 
gauge invariance leading to a calculable, renormalizable theory.' 


However, because of the presence of arbitrary Yukawa couplings to 
calculate fermionic masses and mixing angles as well as arbitrary 
values of coupling constants, vacuum angles and CP violating phases, 
the need for the number of free parameters is excessive (about 18).” 
Also the asymmetry created by the chiral structure of the weak 
interactions and the contrived Higgs sector to break the electroweak 
symmetry has prodded most theorists to believe that there has to be 
a more fundamental theory. Alternates to the standard model are 
_G.U.T. unification leading to large fermion multiplets and fermion 
mass relations,’ technicolor theory suggesting composite structure 
to the Higgs sector,* Supergravity embodied in SU(5) supergravity 
or SU(5) X U(1) supergravity’ reducing the number of free 
parameters to 5 and 3 respectively, and superstrings with only a 
minimum number of free parameters leading to unification of gravity 
and particle interactions at or near the Planck scale.” The other 
attempt at unification involves the notion of compositeness 
suggesting that in analogy with atoms, nuclei and hadrons that all 
..particles (quarks, leptons, gauge bosons and Higgs particles) have a 
finer layer of structure.’ We also believe that at some scale the 
arena of space and time becomes discrete and grainy-like in 
structure. When a discrete quantum theory is imposed on a 
composite system, it turns out that certain measureables reveal the 
composite structure that otherwise would lie hidden beneath a 
quantum formalism embedded in the continuum. ‘This is the 
primary purpose of the study, not to discuss the rigor or “lack of it" 
portrayed by the particular discrete quantum theory that we use,* ° 
but rather to emphasize that independent of the formalism certain 
revealing features of eomposite particles emerge from a discrete 
quantum theory. 
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2. Discrete Time Quantum Theory and Elementary 

Particle Substructure 

The motivation to study the composite structure of quarks, 
leptons, gauge bosons and Higgs particles emerges from the 
enormous success of SU(3) (Flavor) to describe the Hadron spectrum 
and the well known fact that all previous systems (atoms, nuclei, 
hadrons) reveal a composite structure.” '' After the initial 
success of the standard model numerous composite schemes were 
proposed usually with a hypercolor binding mechanism to keep the 
preons together.”  ‘* © 7 Also since quarks and leptons 
are light relative to the composite scale (A ~ 100 Gev) there must be 
a "chiral flavor" symmetry or supersymmetry protecting the quarks 
and leptons from acquiring large masses."* Two of the most 
popular composite models are the Rishon Model of Harari’? and the 
scalar-fermion model by Fritzsch and Mandelbaum.” The Rishon 
model builds quarks and leptons out of three preon composites with 


only two fundamental preons T, V, (7{¢ = 3} Viq = 0) The 


Fermion-scalar model considers a fermion doublet 


two scalar bosons, one of charge , the other of charge -5 ; 


The doublet coupled to the scalar boson of q = generating 


2 
6 


quarks, the doublet coupled to the scalar of q = - generates the 


Es 
2 
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leptons. Higher generations add additional hypergluons. For the 


Rishon model higher generations are built by adding (3T +37) to 


each successive generation. It is clear from these two models also 
that the number of preons is generation dependent. To imprecise 
probes it would be difficult to probe for the number of preons, 
however if the internal dynamics of the preons (constituents of 
quarks, leptons and gauge bosons) is a discrete time quantum theory 
. it-turns.out that observables calculated for the composite system will 
in the zeroth approximation yield the usual Q.M. values with 
corrections expressed in terms of a discrete time interval and 
composite parameters. If the discrete time interval is a random 
variable being the left over fragment of pregeometric "quantum 
correlations" then to calculate the discrete time interval for the 
composite system we can apply the central limit theorem expressing 
the composite discrete time interval in terms of a fundamental 
discrete time interval and the number of preons. For a spin 
-precessing lepton”’ we have calculated the precession frequency in 
‘an external magnetic field, we find , 


oe ae (2.1) 
T 


when p=£h8 


2m 
(t = discrete time interval). 
Here the probability of up and down spin is equal at t = 0. 
Eq. (2.1) gives 


23342 

eB eB? 
bee 

m 96m> 


w « tv (2.2) 


(e = charge of lepton, m = mass of lepton, B = z component 
magnetic field). If there are n preons in the lepton we have from 
the central limit theorem 


T= MT), (tT, = fundamental discrete time interval) 


and 


(2.3) 


Eq. (2:3) would: give an n* dependence of the precession frequency. 
Thus measuring w for different leptons would ascertain the number 
of preons in each lepton. We have also applied this idea to the decay 
of 


(K,, K,), (Bs, B,)-- -(X,, X5) 


using a discrete time asymmetric quantum theory,” the result for 
the short. and long level states is 


8. Ht “is - He s| | 
a - sin if ———- | tan -— |et 
= 1 2h 2 ’ 
Wis = | e * € £ (2.4) 
L,S 


here 


t = discrete time interval, 
time asymmetry parameter. 


mm 
tl 


Both t, € should depend on the number of preons in the 
system (for instance $d, Sd ) and thus very precise 


‘measurements of the lifetime should reveal 


T= Nt, € = Ne, ort = ¥nty, € = YNE,, 


The relations t+ = Vito, € = yne, follow if t and € represent 
the width of a random distribution corresponding to n random 
variables. The relations t = Nt), € = Ne, followife, t are the 


average over a random distribution of n random variables. 

In addition to the above probes to discreteness and composite 
structure we have applied discrete time spin polarization precession 
to a spin one gauge boson with internal (2 preon composite) 
structure. The result was a combination of two sinusoidal functions 
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for the <S,> with the frequencies depending on the internal 
coupling parameters.” Also for elementary particles with internal 
dyon-fermion properties” and elementary particles with heavy 
preons moving in extra dimensions” within the composite particle 
we have calculated spin flip frequencies in an external magnetic field 
and demonstrated that they are sensitive to both the discrete time 
interval t and the internal composite structure mentioned above. 
The models discussed in Ref. (21, 22, 23, 24, 25) all represent 
possible. probes to both the discreteness of time and composite 
structure, though particular simplified systems were used the 
general ideas have far ranging applications. 


3. Markov Environmental Effects on a Quantum System 
In the famous paper of Einstein, Padolsky and Rosen” it was 
emphasized that quantum theory is incomplete and elements of 
reality were missing from the theory. Quantum theory is based on 
"subjective probability" while classical theory is based on objective 
determinism. Whatever interpretation of Q.M. one subscribes to ™ 
“28, 29. 39 we must reckon with the fact that whenever probability 
appears in the interpretation of the theory there may be a composite 
probability theory beneath the observable world or perhaps Markov 
processes adjoined to the normal laws of probability. In this regard 
we have assumed that the Schrodinger theory on the preon or 
substructure level is adjoined to Markov processes that represent 
environmental effects possibly generated by fluctuations of the space- 
time continuum related to the original pre-geometric origin of space- 
time." “ * As we will express below, Markov processes 
represent the correction to a continum Q.M. that a discrete non- 
causal world generate at the present stage of cosmic evolution. 
Perhaps in the very early primitive stages of cosmological evolution 
they were the dominant quantum effects operative before time 
became realized as a continuous variable. It might also be that the 
time scale set by an individual Markov jump might be different than 
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that of the usual time and thus two independent time-like variables 
would appear in the theory. As an example consider a spin /2 system 


in a z component magnetic field with energy E, = + 5 8B ,(q 
m 
= -e) 
The solution to the Schrodinger equation is 
iE,t GE t 


bd -C,ae *+C,Ppe * oo 


Here |C,|?+|C,[? = 1. If an addition to the Schrodinger 


dynamics we have a Markov process operative as described by the 


transition matrix 
. + : —, ~ 7 - 
| l-q q | (3.2) 
p I1-p 


then after n steps if the initial probability of up is ’2 and down is ~ 
Yo we have™ 


~~ 

oo 
+ 

~ 
1) 


(1 
7 Po ett-p-ar (5-2 | 
p+q 2 pt 6%) 


= 
! 

x“ 
I 


- 1 -p-ay(5-— | 
pt+q 2 ptq 
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Here 

P(+), = probability of spin + ¥2 after n steps, 

P(-), = probability of spin - 2 after n steps. 

Now in Eq. (3.1) we replace C, = yPC +),, C, = fP(-),; 
we have 


iE,t (Et 


v= (Phe has ue (3.4) 


giving 


<S,> = h/P(+), P(-), cos af (3.5) 
m 


Then if p # q, the amplitude of <S,> has a certain chaotic behavior 
associated with n. Also we expect p, q to depend on B for high 
magnetic fields. For low magnetic fields we have p = q and 
<5,> = cos (£2 It (3.6) 
Z m 


which is the usual expression for <S,>. Thus a chaotic behavior of 
<S,> and a final state amplitude of 


i [P(*),PC-), - »,|—?4 


(p+qy 


would signal Markov environmental effects. If we have a spin 1 
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gauge boson 
(q = -1) made of two preonic fermions, we have (M,, = mass of gauge 


boson) 


Here 
e -@ 
H = 2(s, +8 (2 is | 
5 1 q = ~@, 
M 
where fe | = > m, = aris for the preons. 
Also 
E = —nB 
M, 
E, = 0 (3.8) 
BOS one 


If we choose 
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ey 
P= /P(+),P(+),e * ao 
-% (3.9) 


+(VPG), PC), «B+ yPC),PC), Bale * 
iE, 
+ (PC) PC)e * BB 


1 
2 


same p, g, then P(+), is same for preon | and 2) we find 


(whereagain P, (+) = > P, .(-)o = = and both preons have 


<§_ +S >= 2ncos 23 
41 M 
? . (3.10) 


x (P25, PO), POPO,] 


If the preons are not the same we still must have 


H- 15 p+“5 B (3.11) 
= —S, oo ; 
m, d Mm, 
€ é . . 
1 2 
eee! ae eee ( = -@€ ’ q — -€,) 
where m, m, 4 ne = but e, ¥ €,,m, # m,. 
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—) = 2 is to insure that the 5. = 0 state has E = 0 
m, m, 


; : e; €, € 
In this case we still have 


(3.12) 


however, the Markov process for preon 1 might have a different p, 


q than for preon 2 because it is seen different by B (external field). 
Thus we write 


iE, 
ee 
Yy = VP,(+), PC ),é » aa 


+(/P,G),P,(- 1),08 +/P,(-), PG), Bale -_ (3.13) 


iE 


+/P,), Pe * BB 


att = 0 
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1 1 € 
Pi( +) = Pal *) = a Pe) BC )p.= 5 (3.14) 
and 
<S>5 = P°(5,, +5.) ? = t 
however at n # 0, t # 0 we have 
S057 a= F (5, +5.) (3.15) 


giving upon using Eq. (3.13) 


JPOEPLCO), Py), | 
+YP,(4), Py (-), P(e 


<§ +S$ >=h cont? (3.16) 
x % M 


+P aPy(n Pa( Dn . 


+ [P,(-)eP,(+), Py), | 


here P,(+),, P,(-), are given by Eq. (3.3) with p = p;, q = q,, and 
P,(+),, P,(-), are given by Eq. (3.3) with p = p,, q = qp. 


In Eq. (3.16) if p, = P2 = 4 = 4 


<5, +9 >= arr (3.17) 


which is the usual expression for a spin | gauge boson. In general 
P.(+1),, P,(-),, P,(+),, Po(-), will vary with n, and B (external field). 
For low fields we expect the expression (Kq. (3.17)) to hold. For 
high fields, however, p,, q;, Pz Q,-Will depend on B. The 
’ experimental problem is to compare Eq. (3.16) with experimental 


measurements of pa ig for discrete steps (n). If the 


expression (Eq. (3.16)) fits the data better than Eq. (3.10) it would 
suggest a preonicsubstructure with different preons. If the formula 
(Eq. (3.10)) fits the data, it would suggest a preonic substructure 
with identical preons. 


The above discussion has outlined probes to both.the discrete 
nature of time induced by pre-geometric fluctuations from Minkowski 
space-time and fundamental Markov environmental corrections to 
Q.M. through discrete time quantum jumps specified by n. If we 
also set t = nt, in Eq. (3.16) then we relate the Markov time 
ordering process to the usual time order specified by t. In two 
previous studies” “ we have outlined a variety of discrete time 
phenomena and Markov environmental processes and in the 
following two tables we summarize the predictions for each model 


for the average spin polarization with <S,>,.,.9 7 ~ for 


fermions and 


<S.>,-,-0 7 4 for bosons.. 


Type of Theory 


Normal quantum theory with no 
Markov influence 


Discrete time difference Q.M. 
with no Markov influence 


Normal quantum theory with 
Markov environmental influence 


Discrete time difference Q.M. 
with Markov environmental 
influence 


md ee 


Table I — 


<S> 


<5,> 


<S,> 


<S.> 


cos| — 
t 


<S,> for spin % fermions (q = -e) 


i] 


TT 


(P(+),, P(.), from Eq, 3.3) 


h eB 


— cos —t (Ref. 35) 


(Ref. 35) 
nY/P(+),P(-), cos( Se 


(Ref. 35) 


hyP(+),P(-), 


(Present Study) 
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Pure discrete time theory with no 
Markov environmental influence 


eBt 


aU ) cos = eFC“ 


<S,> => 


(1-$) +(- 1)" cos = 
(> = €, +€,) 


Ref. 35) 
nt 
ae ae i » €, € are asymmetry parameters necessary 


to specify initial time solution) 


Table II 
Type of Theory 
Normal Q.M. behavior with no <S.> for spin 1 gauge boson 
Markov environmental effects (q = -e) (P(+),, PC), from Eq. 3.3) 
on 2 preon composite gauge 
boson (identical preons) <S, ee eB, (Ref. 36) 
Normal Q.M. with Markov M 7 


environmental effects on 2 
preon composite gauge boson 
(identical preons) 


<S,> = 2heos <* eB, 


[yPCRPCO,+(POEPO,) 


(Ref. 36) 
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Normal Q.M. plus 
environmental effects on spin 1 <S> = h cos eB t 
x 


gauge boson with no composite v2 M, 
structure 


x [2,/P(+), P(0), +2 /P(-), PO), | 


P(+),, P(O),, P(-), = probability of +, 0 -, states calculated 
from Markov chain with initial probability of 4 %, % 
respectively (Ref. 36) 


Discrete time difference Q.M. 

with no Markov environmental <S.> = h cos(a, ~a,)t 
effects (2 preon composite 2 

structure of gauge boson) 


+ eos(a,~a,)t 


(a;, &, ag = parameters depending on internal state of two 
preon composite with S, =. 1, 0, -1) respectively) 


Ref. 36 
Discrete time difference Q.M. / +) P(- = 
with Markov environmental PC Jn ( ),eos(a, a, )t 
effects (2 preon composite <5,> = 2h 
structure of gauge boson with ee [P( on y P( +) cos (a, -a,)t 
identical preons) ; . ; 


Ref. 36 
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6. Normal Q.M. with Markov 


: z om 
environmental Bicct on 2 PCs), P,( + ). P,( =) 
preon composite, ("different 
preons’) See ee ee ee 
2 
+P,(+),P\(-),P,(+)° B 
<S$> =h x cos! 1 


5 M 
+P,(-), Pi(+),P,(-), ? 


: 5 oe: 
+¥P,(-),P2(+),Po(-), 
Here (P,(+),, P,(-),, Po(+),, P,(-), are different for two 
different preons 
Present Study 


In Table I, row 2 would be a probe to both discrete time 
difference effects and the composite structure of leptons if tT = nty. 
‘Row 3 in Table 1 would be a probe to Markov environmental effects 
on normal quantum mechanics. Row 4 would test for both discrete 
time difference effects and Markov environmental influences which 
may or may not include the relationship Tt). = Tarkoy (or the 
interval in discrete time difference effects may or may not be equal 
to the interval for Markov jumps). Row 5 in Table I would test for 
a pure discrete time (not discrete time difference theory) theory and 
predicts chaotic fluctuations in <S,> as evidenced by the term 
(-1)" (e, + €,). 

Table II for row 2 and 3 would provide a test for the compositenss 
of a spin 1 gauge boson since the theory predicts distinctly different 
variations of <S,> with n. Row 4 would provide a test for both 
discrete time difference effects and the composite structure of gauge 
bosons since a,, a, a, are sensitive to the internal preon parameters. 
Row 6 would to provide a probe to test for Markov influences and 
the compositeness of gauge bosons with 2 different types of preons 
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comprising the substructure. 


4. New Frontiers for Discrete Time Q.M. 

As the discussion above has suggested, discrete time notions 
emerge from either pre-geometric theories of how the continuum 
evolved from a discrete set’ or from a fundamental quantum 
uncertainty principle applied to space and time themselves.* 
Probing the discrete property of time might involve very high 
external fields or very rapid physical processes where deviations from 
the continuum may be seen. In addition to spin polarization 
precession, elementary particle decays” that have an "energy 
dependent width-energy relation" generated by discrete time effects 
may also be a probe to discreteness. In an unpublished note” I 
have emphasized that spin correlations of particles emitted in E.P.R. 
type (two particle correlation experiments) may also provide a 
sensitive probe to discrete time effects. Whenever short time 
_ intervals appear in physics whether in particle decays, resonance 
‘phenomena, tunnelling times etc. unaccountable corrections may be 
ascribed to Markov jump processes and discrete time phenomena in 
general. It is not in fashion at the present moment to suggest 
experiments that would probe the discrete nature of time simply 
because there are so many other uncertainties in Q.M. and particle 
theory in general. The existence of higher rank gauge groups, of 
supersymmetry and coupling constant variations all would generate 
corrections to particle life-times and decay asymmetries that 
experimental deviations from the Standard model would be more apt 
to be identified with. This is because of the present "disposition" 
enjoyed by contemporary physics. However, the history of physics 
is filled with examples where a new and stunning discovery came at 
a time when popular belief seemed to discourage it. Such examples 
include the discovery of special relativity, the breakdown of parity 
conservation and the explanation of the Planck radiation law in 
terms of new and radical quantum ideas. In Einstein’s famous book 
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he stresses’ "that the continuum has legitimacy in that it "serves 
to represent the complex of our experience", he also points out that 
any absolute attachment to the notion of continuum physics has no 
justification since it is based on notions outside of our experience. 
It is known that special relativity is violated by discrete physics,” 
but gauge invariance is not. Gauge theories might be the accidental 
result of physics observed in the range (10 TeV > E) emerging from 
the local invariance under a gauge group of a multiplet of fermions. 
‘D, Park has reminded me of the fact that Feynman diagrams and 
even “gauge invariance” are in a certain sense only valid in 
perturbation theory.“ What about the case when photons and 
electrons are inseparable, or quarks and gluons are inseparable, 
doesn’t something deeper have to replace Q.E.D., Q.C.D. and U(1) 
and SU(3), invariance. Perhaps in this "as yet” al domain 
new concepts of space and time will emerge. 

I'd also like to point out that ve the ‘hace of the 
"Thermodynamics of Irreversible Processes” a discrete correlation 
time is discussed that in some sense:represents the "time scale" for 
the principle of microscopic reversibility to hold. This correlation 
time must be bigger than the characteristic molecular interaction 
time and less than the characteristic time for the decay of 
fluctuations. Using this analogy, perhaps the discrete time interval 
in quantum theory is a measure of the time scale for interaction 
between the particle and the background, being greater than the 
preon-preon interaction time and less than the decay time for 
fluctuations from the Minkowski’background. Quantum random 
walks,” “ Quantum Theory on graphs,” the quantization of 
discrete classical systems“ and the search for a discrete graphical 
picture of pregeometry” * all represent pioneering attempts to 
model the world from a discrete set of primitive objects. Although 
there has been a sentiment against modelling Q.M. by Markov 
processes” I feel the systems chosen are to naive to draw any final 
conclusions regarding the underlying probability laws governing 


4r on 
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Quantum theory. To this day, It think it is fair to say that quantum 
theory is still a theory of "subjective probability", the idea of a 
pregeometric world with quantum correlations more fundamental 
than space and time offers” us both a more fundamental picture of 
space and time along with an explanation of how fluctuations from 
the "present" continuum can be interpreted in terms of discrete time 
effects. Let us with a spirit of faith look for instances where 
unexplained corrections to standard physics may be explained by 
discrete physics, thus revealing a new mathematical and physical 
picture of the world. 
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Abstract 

The inverse problem is ”the second half” of quantum mechanics which 
is not widely known yet. The pictures illustrating the fundamental inter- 
relations of spectral and scattering parameters with the details of potential 
perturbations are presented here. We unify the technique of the exactly 
solvable models, program INTERQUANTA that proved itself very useful in 
the picture demonstrations of the direct problem, and the original qualitative 
theory explaining quantum systems by splitting the observables and potential 
perturbations in elementary and universal constituents (atoms, bricks”). 
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INTRODUCTION 

The approach of the inverse problem theory and supersymmetrical quan- 
tum mechanics provides potential perturbations that cause changes of any 
chosen spectral parameters (SP) (without altering other ones). The new 
language [1-4] that allows universally and concisely explains these potential 
perturbations opens the way to new success in theory and applications. 

The new pictures of the shifts of arbitrary energy levels on the energy 
scale and corresponding transformations of potentials and wave functions 
. are given. here.. It is done for different initial potential wells (of finite and 
~ infinite depth, of different shapes): rectangular, oscillator, and multi-step- 
type. Particularly, we shall present the amazing phenomenon of effective 
annihilation of states by degeneration of neighboring energy levels [5]. 

The reader will get a clear notion about the significance of the funda- 
mental spectral parameters (weight or norming factors) which together with 
energy levels form the complete set for one-dimensional systems with the 
purely discrete spectrum. The physical sense of these parameters is that 
each of them relating to a proper bound state serves as a lever controlling 
the region of the main spatial localization of the corresponding state. 

In addition to the already published [1-4] pictures illustrating the con- 
struction of quantum systems with the desired properties it has sense to 
show here more perfect and instructive figures that were derived by the IN- 
TERQUANTA technique which was elaborated by the creation of picture 
books of quantum mechanics [6]. We shall give also explanations of every 
_ partial flexure (sages and arches) of the corresponding potential perturba- 
tions and mention some open problems. 

SHIFTS OF LEVELS IN SQUARE WELLS (infinite and finite), 

There are simple rules of spectral management. The potential pertur- 
bation which shifts only the chosen energy level up (or down) must have 
repulsive hills (attractive wells) in the vicinity of local maxima modulo of 
the corresponding wave function (where the eigenstate is the most sensitive 
to perturbations). And there must be compensating potential wells (hills) 
near the knots of the wave function which keep all other energy levels at the 
previous positions. The main corresponding formulae are given in Appendix 


a 
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Fig.la. Gradual shifts up of only the second energy level: AF, =t=0, 1, 2, 3 
in the infinite rectangular well. Transformations of the four lower bound states 
and of the flat bottom of the initial potential are shown. 
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Fig.lb. The same as in Fig.la, but for t = 4, 4.9, 4.99, 4.999 when Eo, B3 
become almost degenerated. 
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In Fig.1 the gradual deformations of the infinite rectangular potential well 
of the width m and corresponding eigenfunctions are shown in the process of 
approaching of the second energy level E, to the third one. The parameter 
of the energy shift AF, =¢t is gradually changed from ¢ = 0 (unperturbed 
system) to t = 4.999 (the case of almost fulfilled degeneration of the energy 
levels Ey — E>). 

The shape of the potential perturbation has two repulsive hills. At the 
beginning these hills push the chosen second state near the maxima modulo 
of ~2(z) where this state is the most sensitive to potential perturbations. 
‘The three wells of the potential perturbation near the knots of w(x) are 
necessary for the compensation of the influence of the potential hills on the 
other energy levels that must remain at the previous positions. It is possible 
to shift only one chosen energy level (keeping an infinite number of all other 
at the same places !) because all eigenstates are orthogonal to one another. 
The states above the degenerating pair ”lose” gradually the half-oscillations 
from both sides as can be seen for the fourth state in Fig.1b with t=4.9; 4.99; 
4.999. When the neighbor energy levels become very close, the degenerating 
states must simultaneously be orthogonal and similar to one another as the 
solutions of the same equation at almost the same energy. These contradic- 
tory requirements can be satisfied if both states disappear in the middle part 
of the system. The wave functions are divided into two parts being gradually 
pressed into the infinite potential walls of the original well. These two parts 
look more and more like one another by modulus for both the degenerat- 
ing states as the solutions of the same equation at approximately the same 
energy. 

If a wave packet is as a sum of the degenerating states (second and third) 
it will be concentrated at the left wall where (x) and #3(x) have the same 
sign and disappear at the right wall where they cancel one another. Taking 
into account the time. dependent phase factors exp(—7£2,37') we can see that 
there accumulates a phase discrepancy between the two states during some 
time. Therefore the neighbor states change their relative sign and the wave 
packet disappears at the left and concentrates at the right wall. And so the 
packet will” jump” from one wall to the opposite one through the barriers and 
the middle well (without appearing even in the classically allowed region). 
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Fig.2. Gradual shifts down of only the second energy level: AE, = ¢ 
—1, -—2, —2.9, —2.99 in the infinite rectangular well. 
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Fig.3. Gradual shifts up of only the ground state energy level: AF, =%¢ = 
1, 2, 2.9, 2.99 in the infinite rectangular well. 
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3, 6, 6.9, 6.99 in the infinite rectangular well. 
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Transformations of the infinite rectangular potential well of the width 7 
and corresponding eigenfunctions by shifting EF, — F, is shown in Fig.2. 
There are two wells to pull the second state down where it has maximal 
probability to find the particle (wave) and auxiliary repulsion in the middle 
to keep other levels at their previous states. The repulsion in the middle 
stops increasing when w(x) and #3(x) become negligible there. So, the 
middle repulsive hill becomes two-humped. 

Fig.3 corresponds to the shift E, — E2 with the repulsive potential per- 
turbation in the middle region and compensating wells near the vertical po- 
tential walls where the ground state is small. Again the middle repulsion 
becomes gradually two-humped. 

Pay attention to the fact that the limiting shapes of the perturbation in 
the middle are the same in both Figs.2,3. After the effective annihilation of 
£2 and Fy all other states form a new complete set of eigenstates. And the 
place where degeneration of the disappeared levels occurs (below 3) is of no 
importance for other states — we get the same middle well. Compare its form 
with the middle well in Fig.1b. The bottom of the well for the remaining 
upper states must be higher. The states above the degenerating pair ” lose” 

_ gradually. half-oscillations from both sides as can be seen for the third and 
fourth states in the Figs. with t=2.9; 2.99. 

In Fig.4 the transformations of the infinite rectangular potential well of 
the width a and corresponding eigenfunctions while shifting L3 — Ey, are 
shown. Gradually with decreasing 3(r) and 4,(x) in the middle region the 
central potential spike cease to rise (becomes relatively lower: there is no 
need to push the state stronger where it becomes not sensitive). The same 
qualitative rule requires three repulsive potential hills to push up only the 
third energy level (and four auxiliary wells to keep all other levels at their 
previous places). The second potential peak stops to grow when the 3rd 

.and the 4th states become small in the middle. The shapes of the limiting 
potentials for the cases E3 — E, and FE, — E3 must be the same except the 
singular behavior at boundary points z = 0 and z = m. There are confined 
degenerate states in infinite narrow wells behind the infinite high barriers. 

In Fig.5 the transformations of the infinite rectangular potential well of 
the width m and corresponding eigenfunctions in the process of shifting Ey — 
Es are shown. Explanation is qualitatively analogous to that of Figs above. 
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In Fig.5 there are four potential perturbation hills which push the fourth 
energy level up and five compensating auxiliary wells. The two potential 
peaks in the middle stop to grow when the 4th and the 5th states become 
small in the middle. As has been mentioned above, the shapes of the limiting 
potentials for the cases LE, — Es; and E; + E, must also be the same except 
the slightly different singular behavior at boundary points (to provide the 
finite shift of the degenerated energy levels). 

The qualitative difference of the shapes of the limiting potentials in Figs 
1-5 can be explained according to the universal algorithms of symmetrical 
: shifting of states in the new complete spectra after annihilation of the de- 
generated states. 

For finite potential wells ‘he wave motion in the bound states is not 
restricted to the finite interval. Fig.6 shows: the deformation of the finite 
rectangular potential well of the width 6 and depth Vo = —6 and correspond- 
ing bound state eigenfunctions in the process of approaching of the second 
energy level EK, to the third one without changing the reflection coefficients 
of the scattering states. The degenerating states have here the possibility 
to separate exactly from one another by going to too in the opposite di- 
rections. .They are carried away by the auxiliary reflectionless soliton-like 
wells (in order to prevent changing the reflection features of all the states 
of the continuum spectrum when they are moving away). Also the peaks 
appearing at the edges of the initial rectangular well preserve the scattering 
(reflective) features of the initial potential steps which have been smoothed 
by the degeneration. In Figs 7 and 8 the deformation of the finite rectangular 
potential well of the width 6 and depth Vo = —6 and corresponding bound 
state eigenfunctions are shown for Fk, — Fz; FE, > E3 ; E3 > E4; Ey > Es; 
Es; — 0; Fy going down without changing the reflection coefficients. The 
value ¢ = 1 for n = 5 corresponds to the upper bound state lifted to the 
continuous spectrum. The value t = —1;n = 1 corresponds to the ground 
state at the bottom of the initial well. 

Pay attention to the differes:1 shapes of these peaks in Figs 6,7 and Fig.8. 
Compare the transformations of potentials inside the interval [0,1] for the 
infinite and finite initial wells. 
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Fig.8. Symmetrical transformations of square potential having finite depth and 
of corresponding wave functions of all five bound states. when only the ground 
state energy level is shifted down and all other states remain with the same spectral 
and scattering parameters. a) EF, — 0; b) FE; + —4H,; c) By + —14E;. Pay 
attention to the difference in shapes of the peaks at the ends of the initial well 
here and in Figs 6 and 7. 
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of states Ey, Ey; Eo, E3; E3, Eq; E4, Es; Es, Eg; Eg, Ez in 


Fig.11. Annihilation 
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Fig.13. Annihilation of states Ey, 2; 
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Fig.9 shows the transformations of the infinite oscillator well and corre- 
sponding bound state eigenfunctions in the process £3 — E,. All other states 
remain at their previous places. Pay attention to the auxiliary carrier wells 
which move in the opposite directions. In the limit they have infinite depth 
and zero width as the delta-wells. And the corresponding degenerating states 
disappear from the main middle well. Fig.10 shows the transformations of 
the infinite oscillator well and corresponding bound state eigenfunctions in 
the process £, + E2; and FE, — E3. The limiting cases of effective annihi- 
lation of pairs of degenerated states Ey, B2; Ey, Es; E3, E4; E4, Es; Es, Ee; 

..Eg, Ez in- the oscillator well are shown in Fig.11. The explanation of forms 
of the resulting potentials is analogous to the case of the effective deleting of 
energy levels by making the corresponding spectral weight factors infinite or 
zero, as it will be shown below. 

Fig.12 shows the symmetrical deformations of the finite multi-step poten- 
tial well and corresponding bound state eigenfunctions in process of approach 
of the third energy level 3 to the fourth one. The unperturbed initial sys- 
tem is‘chosen with the equidistant discrete spectrum. It is interesting that 
the equidistancy requires the initial multi-step well to be approximately of 
oscillator shape in the range of the lowest part of the spectrum where their 

‘energy. levels coincide. Compare the shapes of the potential perturbations 
(find the common and different features) in Fig.12 and Fig.9. Pay attention 
to the potential peaks corresponding to every step of the initial well (for 
conservation of the reflection in scattering states by potential deformation). 

Let us mention two kinds of the universal form of the auxiliary wells which 
carry the degenerating bound states out of the main well (for transformations 
of infinite and finite wells). Every carrier well takes away only half oscillation 
of the initial wave function. The oscillations which remain in the main well 
disappear in the limit of degeneration. 

In Fig.13 are shown the deformation of the finite multi-step potential well 
-and corresponding bound state eigenfunctions in the limiting cases: BK, — Eo; 
E, — E3; Ey — E3; E3 — E4; Eq — Es; Es — Ee; Eg — Ez. Compare the 
shapes of the perturbed potentials in Figs 13 and 11. 

SHIFTING BOUND STATES IN SPACE (Changing spectral 
weight factors) 

The spatial shifts of any chosen states are performed by auxiliary wells 
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. Fig.14. Shifting of the ground state from the oscillator well. 
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Fig.15. Shifting of the 2d and 3rd states from the oscillator well. 
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of the perturbation potentials carrying the states to be moved and analo- 
gous to that considered above (for shifts of energy levels). This confirms 
the economy of the algorithmical tools of the qualitative theory of quantum 
spectral control. There exist favorable conditions for retaining the standing 
wave in the carrier wells (for a half oscillation only) for the chosen energy. 
The corresponding formulae are given in Appendix 2. 

Fig.14 shows the transformations of the infinite oscillator potential well 
and the corresponding bound state eigenfunctions in the process of increas- 
ing the spectral weight parameter (pre-exponential factor in the asymptotic 

‘behavior) M;, of the ground state.: Ground state is moved away from the 
main well by a special auxiliary carrier well. There exist favorable conditions 
of retaining the standing wave in the carrier wells (for a half oscillation only) 
for the chosen energy. All other states have their last half-oscillations inside 
the carrier well but with the disappearing amplitudes. 

Fig.15 shows the gradual removing of the second and third states. Only 
one half-oscillation of the standing wave is carried away. The other oscilla- 
tions remain in the main well but they have negligible amplitudes (unfavor- 
able conditions for the standing wave there), so that the whole normalization 
of the removed state is concentrated in the carrier well. All other states re- 
main in the deformed initial well and have their last half-oscillations inside 
the carrier well but with the disappearing amplitudes (destructive interfer- 
ence). The main well in the limit becomes narrower in comparison with the 
initial one above the removed state because every upper state lost a half- 
oscillation and must be shorter. The additional deepenings at the bottom of 
the main well are needed to keep at the same position the energy levels of 
the states below the removed one (to pool them down at the previous places 
when the narrowing of the main well pushes them up according to the rule 
considered above). 

DELETING OF THE CHOSEN ENERGY LEVELS 

In the limiting case M,, — oo we have the effective disappearance of the 
nth energy level from the bound state spectrum. 
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, Fig.16. Deleting of states Ey; Bo; E3; E4; Es; Ee in oscillator well. 
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. Fig.17. Deleting of states E42; E24; E49; Egg in oscillator well. 
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Fig.16 shows the limiting cases: deletion of states Fy; 2; F3; E4; Es; Eg. 
As in Fig.15, all the states above the removed one lost one half-oscillation 
and became shorter. To keep at the same positions the energy levels below 
the removed one, we need to pool them down with the proper number of 
partial wells at the bottom. The most sensitive is the nearest neighbor of 
the removed state. For more deeper states the narrowing of the upper part 
of the main well is not so important. So the shape of the bottom is mainly 
determined by the forces applied to this state. 

Fig.17 shows the limiting cases: deletion of states E12; Fo4; E49; E99. Pay 
attention that the frequent potential oscillations do not alter significantly the 
wave functions of the lower states (approximately zero average influence). 
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APPENDIX I 

Symmetrical deformations of initial potential wells and bound states by 
shifting only the chosen nth energy level are given below [7]. 


v(x) =0 (x) — 2{Inw(2)}", (1) 


where 
w(2) = [v(x),b (2, En)]. (2) 
In this expression symbol [ . , . ] denotes Wronskian: 
[w(x),b (2, Bn)] = w(x) '(2, En) — w'(x) P (2, En), (3) 


where 7 (x, #,) is the n-th (normalized) bound state wave function of the 
reference system. As w(x) one should take the solutions for initial symmet- 
rical potential such that their parity must be opposite to one of the state to 
be shifted. For example, in the case of oscillator v(x) = x? hypergeometrical 
analogues of sine and cosine should be used: 


w(x) = exp(—2?/2) iFi(1/4 — (Ex + t)/4,1/2,27), n is odd, 
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w(x) = exp(—27/2) 1.F4(3/4 — (En + t)/4,3/2, 27) a, n is even. 


And for the rectangular well as w(x) should be taken either sine or cosine 
(in the case of the well of finite depth these solutions, of course, must be 
smoothly sewed together with exponentially decreasing solutions outside the 
well). The expression for eigenfunction at the shifted energy E, + t is given 
by 

(x, Ey +t) =p (2, En)/w(z). (4) 
In the case of the whole line eigenfunctions for other (unshifted) energy levels 
Em are written as following (% 4 m) : 


—1(1/w(z))w() | (y, Bn) B (y; Bm) dy. (5) 


{ 
8 — 5 
a 


APPENDIX II 
Spectral parameters which determine the boundary (asymptotic) behav- 
ior of bound state wave functions are the levers of control of the main space 


localization of the chosen states. Let 7 (x) be an eigenfunction of the refer- 


ence Hamiltonian 0 (x) = 2’ related to wth energy eigenvalue, M and M be 
the factors by the decreasing exponential asymptotical tails ~ x” exp(—zx?/2) 


of (xz) and ib (x) respectively. If the problem is to alter only the weight fac- 


tor M of chosen state provided that all other spectral parameters remain 
unchanged, then the transformed potential is 


° d 
V(z) =v (2) — Ty (#12): (6) 
Here expression for K has the form 
(1 — M?/ M*) ¢ (x) ¥ (y) 


K(t;y)= eS 
1—(1— M?/ Mm?) f b*(y)dy 
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Formula for the eigenfunctions of new Hamiltonian at unaltered energies E, 
reads 


a 

= 

a 

Il 

e 

= 

+ 
acs 

= 

Go 


y)-by (y)dy. (8) 


It should be noted here that similar formula for K is valid also for varying 
c in G-L-approach (as well as in the case of finite interval). We must only 
replace c by M and the integration limits in the denominator (from 0 till z). 
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